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ABSTRACT 


Computer analysis with IDL was performed on polarized infrared image data 
obtained during the MAPTIP coastal experiment. Analysis of contrast enhancement 
of the oceanographic ship Hr. Ms. Tydeman and a Lynx helicopter using 
polarization filters is presented. Contrast of vertical and horizontal polarized 
images was determined by two methods. Unpolarized image contrast was 
developed from polarization data and used to calculate contrast improvement 
factors for vertical and horizontal polarizations. Background effects upon contrast 
were explored. Dependence of projected ship area on contrast improvement was 
analyzed. Vertical and horizontal polarized images were compared to verifiy 
preliminary observations of negligible predominant polarization of man-made 
targets. 

Results of this analysis confirm that vertically polarized sea background 
radiance can be reduced by polarization filtering ts enhance ship-to-backround 
contrast. Man-made targets appear to present negligible polarization. Sky 
background immediately above the horizon also appears unpolarized. Combining 
radiance from both sky backround above the horizon and sea background to form 
an average background reference radiance reduces overall contrast improvement. 


Ship-to-background contrast should be calculated with respect to the horizon 


location. 
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I. INTRODUCTION 


As the U. S. Navy moves to increase weapons systems capability in the littoral 
environment, the ability to detect craft with low or negligible radar cross-section such 
as fiberglass hulled gunboats demands increased consideration. Infrared detection and 
tracking systems enhance the combat capability of both ships and aircraft in this 
environment. However, infrared signatures of interest often exist in substantially 
cluttered backgrounds. Design of infrared systems capable of satisfactory operation 
despite clutter requires improvement of target to background contrast to meet the 
system detection threshold limit. Successful clutter resistant infrared systems could 
fill the gap between longer range radar detection and short range daytime visibility 
detection. 

Current infrared detection and tracking systems use a variety of filtering 
methods for target detection and recognition. Various spatial and spectral filtering 
techniques achieve some improvement in target to background contrast. However, 
polarization filtering could provide significant improvement in target to background 
contrast. 

Polarization filters reduce polarized apparent sea radiance clutter in the target 
scene and improve target to background contrast. Specifically, the electric vector of 
sea radiance clutter appears perpendicularly (horizontally) polarized within sun glint 
regions and parallel (vertically) polarized otherwise (Cooper, Crittenden, Milne, Moss, 
Gregoris, 1992). Man-made targets appear much less polarized than sea backgrounds. 
Polarization filters aligned horizontally while viewing target scenes outside of sun glint 
indicate a greater reduction of sea background radiance clutter compared to vertically 
aligned filters viewing the same target scene (Cooper, Lentz, Walker, Chan, April 94). 
For infrared system applications, the target scene contrast achieved with polarizers 
should exceed that of an unpolarized scene in order to reduce the system detection 
threshold value. This in turn would increase minimum detection and recognition 
ranges. 

This thesis analyzes marine environment polarization images taken during the 
MAPTIP experiment of 18-31 October 1993 off the coast of the Netherlands. These 
images provide the apparent radiance from target scenes involving ships and aircraft 


as viewed through vertical and horizontal filters. Computer code was generated to 














perform the analysis, providing the apparent average ship and background radiance in 
vertical polarization, horizontal polarization, and unpolarized representation. Also 
provided are vertical, horizontal, and unpolarized target to background scene contrasts. 
As a result, contrast improvement attained through vertical and horizontal 
polarizations compared with unpolarized scenes is produced. 

Chapter II provides a discussion of essential infrared radiation fundamental 
concepts and equations necessary for more involved development of infrared radiation 
theory. Chapter III describes the concepts involved in light polarization and the theory 
of infrared polarization in the marine environment. Chapter IV describes the image 
data utilized and also the specifics of the analysis procedure. Chapter V provides a 


discussion of the results, conclusions and recommendations. Appendix A contains the 


three computer programs used and Appendix B provides the results of the analysis. 








II. INFRARED RADIATION FUNDAMENTALS 


A. THEORY OF INFRARED RADIATION 


Infrared radiation is a part of the total electro-magnetic spectrum. The 
infrared region spans wavelengths from 0.72 pm to 1000 pm as shown in the diagram 
of the electro-magnetic spectrum in Fig. 2.1 (Grum and Becherer, 1979). The red end 
of visible light borders the infrared region's lower limit and the microwave end of radio 
waves borders the upper limit. The infrared region subdivides further into near (0.72- 
1.5 ym), middle (1.5-5.6 pm), and far (5.6-1000 pm) infrared regions (Hackforth, 1960). 
Any material at a temperature above absolute zero emits radiation. A body above 
absolute zero contains a finite value of internal energy associated with the excitation of 
atoms and molecules. Excitation comprises the storage of energy in electronic, 
vibrational, and/or rotational states. The levels of excitation in vibrational states occur 
in multiples of some discrete value of energy. When the atom or molecule returns to 
the lower energy state, it releases radiation in the form of a photon at a frequency 


proportional to the change in energy between states, 


AE=hv 2.1 


where, AE = change in energy between states (J or eV) 

h = Plank's constant (6.6256 x 10™ J-s) 

v = frequency of released radiation (Hz) 
The change in energy associated with electronic states typically ranges from one to ten 
electron volts (1.24-0.124 pm), whereas vibrational and rotational state changes in 
energy are typically 0.1 eV (12.4 pm) and 10% eV (12,424 pm), respectively. Therefore, 
molecular infrared radiation primarily involves vibrational changes of state with 
rotational changes superimposed. Fig. 2.2 shows the absorption spectrum of diatomic 
molecules of HCL gas. The actual emission spectrum contains rotational changes of 


state superimposed upon the vibrational changes. (Cooper 4253 notes, 1994) 








pra Optical —_——_—| 


infrared 


Im 1pm idem Oicm icm cm im = fm = =Om thm km (0Okm WAVELENGTH 
107 3x10 nad? a0 3x00 3x10 32103 3x02 iio! 3x10 310° 3x00® 3x10” 31106 32105 32004 3x10? — FREQUENCY, Hz 


{Middle infrored} Far infrared 


30 WAVELENGTH, pm 
WAVENUMBER; cm! 





Figure 2.1 Electro-magnetic spectrum (Grum and 


Becherer, 1979). 
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Figure 2.2 Vibration-rotation spectrum of HCL (Cooper 4253 notes, 
1994). 


1. Infrared Definitions and Units 

Several definitions and units require introduction for an understanding of 
infrared science. Table 2.1 shows fundamental radiometric quantities associated with 
infrared radiation. Radiant exitance (emittance) or flux density, M, is the radiant flux 
exiting an infinitesimal area divided by that infinitesimal area. Irradiance or radiant 
flux surface density, E, is the radiant power incident upon a unit area of surface. 
Radiant intensity, I, is the radiant power exiting a point source along a given direction 
within a unit solid angle. Radiance or radiant intensity surface density in a given 
direction, L, is the radiant power per unit solid angle exiting from or incident on an 
area projected perpendicular to the direction of radiant energy flow. (Grum and 
Becherer, 1979, Cooper Class Notes, 1994; Wolfe and Zissis, 1978, Bramson, 1968, 
Smith, 1993) 

2. Blackbody Radiation 

The blackbody is fundamental to the study of infrared radiation. Figure 2.3 
shows typical blackbody radiant emittance curves for several temperatures. The total 
radiant emittance or area under the curve and the peak value of emittance depend 
upon the blackbody temperature, the peak value shifting to the left at greater 
temperatures. The radiation emitted at a given temperature is a maximum for a 
blackbody. As the standard for radiation theory, blackbody emission and absorption 
efficiency, designated as the emissivity factor, equals one.(Hackforth, 1960) 
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a. Planck's Law 
During the late 19th century many scientists grappled with deriving the 


equation fitting the experimentally obtained radiant emittance curves similar to those 
in Figure 2.3. Max Planck successfully established this equation and simultaneously 
ushered in the era of quantum theory. Planck theorized that oscillating modes of 
frequency w excite higher oscillation energy levels in discrete increments or quanta. 
This theory led to Planck's law which follows, 
2nc7h 
he 


M,(T)= [Wem™*cm~"] 2.2 


a5(e ®T_1) 


where M,(T) = spectral radiation emitted into a hemisphere per unit area per unit 
wavelength interval (Wem?cm") 
T = absolute temperature of the body (°K) 
X = wavelength (cm) 
h = Plank's constant (6.6256 x 10“ J-s) 
speed of light in a vacuum ( 2.998 x 10” cm/s) 
Boltzmann's constant (1.38054 x 10” J-K"). 
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As seen from this equation and Figure 2.3, as the temperature increases the radiant 
emittance increases significantly. In fact, the total radiant exitance increases with the 
fourth power of absolute temperature as shown by the Stefan-Boltzmann law described 
below.(Hudson, 1969, Hackforth, 1960) 
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Figure 2.3 Radiant emittance (Hudson, 1969). 





b. Emissivity 
No real bodies emit and absorb radiation ideally as the blackbody. 
Emissivity depends on the body constituents, temperature, incident radiation quantity 


and wavelength, and surface finish. The ratio of a body's total radiant emittance, M, 








to that of a blackbody, Mgz,, defines the body's emissivity, 


€ -M 2.3 
Mp, 


Again, blackbody emissivity equals one. Some materials approximate blackbody 
behavior. Many others, called gray bodies, display emissivity factors less than one, but 
are independent of wavelength as is a blackbody. Still other materials, called selective 
emitters, emit and absorb radiation variably dependent upon wavelength. The 


following equation accounts for the possible spectral dependence of emissivity, e(), 


fequm,aa 
_ 0 
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- 2.4 
[maa 
0 


Figure 2.4 displays radiant emittance curves for a gray body and selective emitter 


compared with a blackbody curve.(Hudson, 1969) 
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Figure 2.4 Radiant emittance for selective, graybody, and 
blackbody emission spectra (Hudson, 1969). 














c. Stefan-Boltzmann Law 
The integration of Planck's Law over all wavelengths produces the total 


radiant exitance emitted from a blackbody as a function of temperature as shown 


below, 
M(D= [ M, (Tdi =oT*[Wicm?} 2.5 
0 
where, M(T) = total radiant emittance at temperature T (Wcm”) 
o = Stefan-Boltzmann constant (5.670 x 107 Wem?K*). 


For a graybody, this function becomes, 
M(T)=€0T*. 2.6 


Quite often, spectral filters used with infrared systems permit only a specified 
wavelength interval of irradiance to reach the detector. Therefore, effective radiant 
emittance then results from integrating M(T) over only the desired wavelength 
interval.(Hudson, 1969, Hackforth, 1960) 

d. Wien's Displacement Law 

Differentiating the radiant emittance M, of Equation 2.2 and setting 
the result equal to zero yields the wavelength of maximum emittance. This is known 


as Wien's displacement law, 


1, TK 2.7 


where, = wavelength of maximum emittance (cm or pm) 


Mn 
T = absolute temperature (°K) 
K = constant = 0.2897 cmK for Am in cm = 2897 pmK for A,, in pm. 

Since the wavelength of maximum emittance times the absolute temperature equals a 
constant, then A,, displaces to shorter wavelengths as temperature increases as shown 
in Figure 2.3 by the dashed line.(Hudson, 1969, Hackforth, 1960) 

e. Kirchoff's Law 

A blackbody absorbs all incident radiation and emits an equal amount of 
radiation. Applying conservation of energy to a non-blackbody yields, 


where U, = incident radiant energy 





U,=U,+U,+U, 2.8 


U, = reflected radiant energy 


U,, = absorbed radiant energy 
U, = transmitted radiant energy. 


Normalizing with respect to U; produces 


1l=pt+a+t 2.9 


where, p = U/JU; = reflectivity 
a = UJU, = absorptivity 
t= U/JU; = transmissivity. 
Kirchoff showed that the ratio of radiant emittance, M, to absorptivity for a body at a 


specified temperature equals a constant value of the radiant emittance of a blackbody, 


M,,, at the same temperature or 


M 
Mp7 2.10 


Substituting in the Stefan-Boltzmann equation produces 


oTt=£9F- 2.11 
a 


a =e, 2.12 


So, the emissivity of a body equals its absorptivity. For opaque bodies, t = 0, thus 


e=l-p 2.13 
which means good emitters reflect poorly.(Hudson, 1969) 

f Lambert's Law 

Since a blackbody is a diffuse radiator, it has the same radiance in all 
directions. However, the radiant intensity of an element of the blackbody surface 


varies with the angle from the normal to the surface. Lambert's Law states this effect 


as, 











dI=LdAcos6 2.14 


where, dI = element of radiant intensity 
L = radiance 
dA = surface area element of body 
@ = angle from normal to observation direction. 
Figure 2.5 shows the geometry involved in Lambert's Law.(Jamieson, Plass, McFee, 
Grube, and Richards, 1963, Hackforth, 1960). The x steradians in Figure 2.5 denotes 
the quadrant of space to the left of the source area and above the horizontal line that 


passes through the center of the source. 


Source $ 
of area A 


Receiver 
R 


Figure 2.5 Lambert's law of cosines (Hackforth, 1960). 





B. INFRARED RADIATION SOURCES 


As stated before, any body above absolute zero emits radiation. Types of 
infrared radiation sources include artificial, natural, and background. Figure 2.6 
depicts the many types of infrared radiation sources. Artificial sources consist of 
primarily man-made materials such as lamps, filaments, and laboratory blackbody 
sources. Trees, water, mud, vehicles, and people makeup terrestrial natural sources. 
Atmospheric gas molecules, dust, and water droplets are also terrestrial sources. 
Terrestrial sources emit, reflect, absorb and reradiate radiation characteristic to each 


particular material. Celestial natural sources include the sun, moon, planets, and 


11 








stars. Terrestrial and celestial sources serve as background source radiation if super- 


imposed on the primary source radiation.(Hackforth, 1960) 
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Figure 2.6 Radiation Sources (Hackforth, 1960). 


components 
® Clouds 
@ Environmental 


Background radiation also differs with respect to day or night. Figure 2.7 


depicts the relative radiant energy of clear day and night skies. During the day, 


reflection and scattering of sunlight by atmospheric particles, water droplets, dust 


particles, and clouds as well as emission from atmospheric constituents at ambient 


temperatures contribute to background radiation. Daytime background radiation 


increases in the direction of the sun and horizon and remains higher than nighttime 


levels in all directions. The angle of elevation and atmospheric temperature affect both 
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Figure 2.7 Spectral energy distribution of background radiation from clear day and night 
skies (Hackforth, 1960). 


daytime and nighttime radiation intensity. At night, the high frequency radiation 
scattered by small particles and molecules vanishes. Background radiation at night is 
emitted from atmospheric molecules and particles with the moon and stars providing 
additional low-level radiation. (Hackforth, 1960) 

Discrimination of a primary source against its background requires a finite 
primary source temperature different from its background, or a difference in source 
and background emissivities. The magnitude of the required temperature difference 
depends upon both the detection system capabilities and the ambient photon noise 
environment. The background signal generates noise for the detector. Systems use 
various methods such as spatial, spectral, or electronic filtering to reduce the 


background signal to enhance primary source or target detection. (Hackforth, 1960) 
C. ATMOSPHERIC PROPAGATION OF INFRARED RADIATION 


Whatever the source, the detectable radiation received at a detector differs from 
that radiation emitted at the source. The atmosphere or some other intermediate 
medium through which the radiation passes attenuates the radiation through 
scattering and absorption. These effects occur independently from the source, but vary 
according to the wavelength of radiation. Current meteorological conditions play an 
important role in radiation propagation. As the atmosphere changes so does its 
propagation characteristics. Cloud type, base, and ceiling, fog, suspended particles, 


and molecular constituents of the atmosphere all impact absorption and scattering 








effects of the atmosphere on infrared radiation. Beer's Law quantifies the total 


attenuation by the atmosphere as 
K(z)=I(O)e 2.15 


where, 1(0) = initial intensity of radiation 
I(z) = intensity of radiation after propagation distance z in the 


atmosphere 
p = extinction coefficient of atmosphere 
Zz = propagation distance. 


Therefore, the transmittance, t results from, 


p=! 2 -, 2.16 


The atmospheric absorption (,) and scattering (u,) coefficients sum to produce the 


total atmospheric extinction coefficient (y) 


B=H,+H, 2.17 
where, 

Bo=*, TK, 2.18 

#,=9,,79, 2.19 


and x,,= molecular absorption coefficient 

kK, = aerosol absorption coefficient 

6,, = molecular scattering coefficient 

6, = aerosol scattering coefficient. 

The magnitudes of each coefficient depend upon the atmospheric composition and 
constituent density and distribution. Aerosol scattering dominates infrared radiation 
extinction below 1.2 ym. Above 1.2 pm, molecular absorption dominates. In addition, 
atmospheric transmittance depends significantly on wavelength.(Cooper class notes, 
1994) Figure 2.8 shows transmittance as affected by infrared absorption by the 
predominant atmospheric molecules over a horizontal path of 6000 meters at sea level 


with 17 mm of precipitable water (Hudson, 1960). 
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1. Absorption Effects 

Absorption occurs when infrared radiation interferes with a free molecule, atom 
or a liquid or solid at about the frequency corresponding to an energy level separation, 
and the electric field of infrared radiation excites that material to the higher energy 
level. The radiation disapears as a result (Hackforth, 1960). Lorentzian curves depict 
the frequency lineshape of absorption which follows the classical resonance 
phenomenon. Inhomogeneous broadening results in frequency lineshape that are 
Gaussian. Also, collisional, Doppler, or pressure broadening increases the possible 
range of frequencies for absorption (Cooper class notes, 1994). Carbon dioxide and 
water dominate infrared molecular absorption in the atmosphere (Hudson, 1969). 

To a much lesser extent, absorption and scattering by aerosols affects infrared 
propagation. Particles suspended in the atmosphere such as dust and volcanic ash 
generate aerosols. Their distributions and number densities vary considerably from 
one environment to another. Several aerosol models available predict aerosol 
composition of the atmosphere and effects on propagation.(Cooper class notes, 1994) 

2. Scattering Effects 

Scattering results from the absorbing molecule, atoms, or aerosol subsequently 
reradiating but over a 4n solid angle, thus reducing the directivity of the original 
radiation. Where the incident radiation wavelength greatly exceeds the size of the 
molecule or aerosol particle, the following reradiation approximates isotopic radiation. 
As particle size increases compared to the incident radiation wavelength, the radiation 
increasingly tends to be scattered into the forward hemishpere. Rayleigh scattering 
theory effectively describes small particle scattering, while Mie theory covers the entire 
range of particle sizes.(Cooper class notes, 1994) 

a. Rayleigh Scattering 

Rayleigh scattering theory describes scattering for small particles of 
radius r, typically molecules, where 2nr/A«1. Other assumptions for Rayleigh 
scattering include nonionized isotropic particles, index of refraction approximately one, 
and incident radiation whose frequency is different from particle resonant frequencies. 
Since diatomic nitrogen and oxygen constitute most of the atmosphere, the isotropic 


assumption requires a molecular anisotropy correction factor. The equation governing 
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Rayleigh scattering follows, 


1 xs? (1 400820) 2.20 


| NS et a 





where, I = intensity of scattered radiation 
I' = intensity of incident radiation 
K = characteristic of medium 
V = particle volume 
r = distance to particle from observation point 
X = wavelength of radiation 


6 = scattering angle. 
Synchronously induced dipole moments from the incident radiation produce Rayleigh 


scattering. (Smith, 1993) The scattered radiation frequency equals the incident 
radiation frequency. Thus Rayleigh scattering is elastic (Coulson, 1988). 

b. Mie Scattering 

Mie scattering applies to both large and small particles compared to 
incident wavelength. This theory assumes isolated spherical particles as the scattering 
medium. Although most particles are non-spherical, large numbers of randomly 
orientated non-spherical particles approximate spherical particle distributions except 
for determining properties of backscatter and polarization. (Smith, 1993) As with 
Rayleigh scattering, Mie scattering describes the excitation of particles by incident 
radiation and the reradiation of that energy into a 4n steradian solid angle (Coulson, 
1988). The Mie scattering cross section per unit volume follows, 


Boea= | Coca (rr 2.21 
0 


where, B.., = scattering cross section per unit volume 


C,,,(r) = scattering cross section 


n(r) = particle size distribution function (# particles/unit volume 
between r and r+dr (Smith, 1993) 
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Figure 2.8 Atmospheric transmittance over a horizontal path of 6000 meters at sea 
level with 17 mm of precipitable water (Hudson, p.115). 














Ill. INFRARED POLARIZATION IN THE MARINE ENVIRONMENT 


A. POLARIZATION OF LIGHT 


Much experimentation by Fresnel and Arago led to the conclusion that a light 
wave consists of two components perpendicular to the direction of propagation. The 


components for a wave propagating in the z-direction may be represented by 


E,(z,t)=E,cos(wt+8,) 3.1 


E(t) =E,,cos(wt+6,) 3.2 


where E,, and E,, denote the component electric field amplitudes, 6, and 6, are the 
respective phases, w is frequency, and t is time. E,(z,t) and E,(z,t) produce a vector in 
the plane perpendicular to propagation that generates a curve in that plane known as 
the polarization ellipse described by 
Ee Ay as teat =sin?d 3.3 
Ey? Ey Fa Fey 


where 6 = 6, + 6,. This ellipse is shown in Fig. 3.1. Degenerate forms of the 


polarization ellipse include circular and linear polarizations.(Collett, 1993) 





Figure 3.1 Polarization Ellipse (Collett, 1993). 








The polarization ellipse is applicable to completely polarized light only, and not 
useful for describing partially or unpolarized light. Sir George Gabriel Stokes resolved 
this problem by describing light with observable intensities. He determined that light 


of any polarization can by described by the four Stokes parameters given by 


S,228,7+58,2+8,7 3.4 
S,=E,,’-Ey, 3.5 
S,=2E,E,c0sd 3.6 
S,=2E,,E, sind 3.7 
where 
Sq Eq, +E,’ 3.8 


(Collett, 1993) Then the beam intensity is given in terms of Stokes parameters as, 


18,6) =515,+S,c0s(20) +S,cos(6)sin(26) +5,sin()sin(26)] 39” 


where @ is an induced phase angle (Collett, 1993). 

Stokes parameters of two independent light beams may be added to obtain 
Stokes parameters for the resultant beam. Thus partially polarized light may be 
represented by the summation of unpolarized light with completely polarized light. 
Both unpolarized light and partially p: :arized light may be represented by summing 
two oppositely polarized beams. For unpolarized light, the intensities of the two beams 
must be equal. Thus a partially polarized target scene in a marine environment may 


be represented by two beams of vertical and horizontal polarization.(Collett, 1993) 


B. REFLECTION AND REFRACTION 


Reflection and refraction of light is described by the Fresnel equations and 


Snell's law. The incident, reflected, and transmitted waves are shown in Figure 3.2. 
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Serving as a starting point are Maxwell's equations of electro-magnetism in a source- 
free medium, 


wai--y* 3.10 
at 
wxli=e" SE 3.11 
at 
ve= 3.12 
wH=0 3.13 


where E and H are the electric and magnetic field intensities, respectively. The 


complex permitivity and complex permeability are denoted by e* and p*, respectively. 
For plane waves of the form, 


E,=E,¢" 


3.14 

i Jat-y 5 
H,=Eye y 3 15 
y=jo(yxe*)*=0 +78 = 


incident at a boundary between two media of different complex refractive indices n, 


and n,’ the angle of reflection equals the angle of incidence. The angle of transmission 
@, as a function of incident angle 0; is given by Snell's law 


n, sin(8,)=n, sin(6)). 3.17 














For nonmagnetic media the Fresnel equations can be written 


__sin(®,-6) 


r 
+" sin(6,+@) 


a tan(6,-8,) 


| tan(6,+6,) 


_ 2sin(6,)cos(8)) 
* — sin(6;+6,) 


___ 2sin(@,)cos(8,) 
| sin(6,+6,)cos(8,~6,). 


The Brewster or polarizing angle occurs where 


"9 
ny 
and 
6,+8,=90. 
The reflection coefficients then become 
n= 
r #0. 


Thus the reflected light is completely polarized normal to the plane of 


incidence.(Sandus, 1965) 
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C. SEA SURFACE INFRARED POLARIZATION 


Background infrared radiation from the sea surface is composed of reflected 
radiation from the sky, sun, moon, stars, etc. Another contribution to background 
radiation comes from thermally emitted radiation from the water, and molecular and 
particle volume scattering (Sidran, 1981). Sea surface reflection is governed by the 
Fresnel equations with polarization of reflected radiation greater near the Brewster 
angle. Polarization of emitted radiation is somewhat more involved. 

1. Sea Surface Reflectance and Emissivity 

Polarized light emanating from a surface may be due to emitted radiation as 
well as reflected radiation. Arago showed that light from an incandescent solid or 
liquid is partially polarized even for non-polished surfaces. Arago's work was 
qualitative. Millikan later described Arago's work in quantitative terms.(Sandus, 


1965) 


Medium 2 


B 
Medium tf 





Figure 3.2 Reflection and Transmission at a Boundary (Sandus, 1965). 








Millikan studied the polarization characteristics of uranium glass fluorescence. 
He observed that uranium glass fluorescence emission reached a maximum 
polarization at near grazing angles. Furthermore, the emitted light was polarized 
parallel to the plane of emission. As described above, at least partial polarization 
should be observed for reflected light. However, this light would be polarized normal 
to the plane of incidence. Therefore, the polarized light Millikan observed from the 
fluorescence of uranium glass could only have come from emitted light undergoing 
refraction.(Sandus 1965) 

Extending the understanding of emission polarization, Basener and McCoyd 
developed a model to predict the percent polarization of emitted infrared radiation 
from the sea surface as a function of wind speed and viewing angle. This model is 
based upon a sea surface of tiny plane facets of varying orientations. The facet 
orientation is statistically distributed depending on wind speed.(Basener and McCoyd, 


1967) 
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Figure 3.3 Reflection geometry for a facet of a wind roughened sea surface (Gregoris, Yu, 


Cooper, Milne, 1992). 
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Given that the sea surface emits polarized radiation, Basener and McCoyd 
investigated the affects on emissivity of the variable nature of the sea surface. Figure 
3.3 shows the reflection geometry for a facet of a wind roughened sea surface. Each 
facet is assumed flat. For a flat blackbody surface z=f(x,y) the directional emissivity 


€,(0,0) is related to the directional reflectance R(6,) through 


€,(8,6)=1-R,(8,6) 3.26 


where r is polarization mode. Thus, the radiant flux emitted per unit solid angle per 


unit area per unit wavelength interval is 


F (0,0)=—,(0,0)e0s(®)B(4,7) 3.27 


where B(A,T) is the Planck blackbody function. The reflectance is found from the 


Fresnel equations above. For complex indices of refraction, the emissivity is given by 


ga 3.28 


" (sz, 
where z , is a complex index of refraction.(Basener and McCoyd, 1967) 

The above equation is for a flat surface and the sea surface is not flat. 
However, the sea surface can be approximated by many facets with distributed 
orientations. For the facets to be considered flat then they must be large compared to 
the wavelength of light. Sea surface waves easily meet this constraint. Referring 
again to Fig 3.3, the vector O points in the direction of the observer, S is the sun 
direction, and the vector N is the facet normal. Therefore, the emissivity depends on 
the angle between S and N. The radiant flux emitted per unit solid angle per unit 


wavelength interval from a unit horizontal area is given by 
1 
F,0,6)= BOD f«,(B,y)sec( BB, y)e0s(#)dea 3.29 


where e'(B,y) is the emissivity referenced to the main plane of vision defined by S and 
k, where k is the unit vector in the z-direction, B is the zenith angle of the facet 
normal, and y is the azimuth of the facet normal. The facet distribution function is 


denoted by f(B,y) and y is the angle between S and N. Cox and Munk (Basener and 








McCoyd, 1967) have determined experimentally that the distribution of slopes given by 
f(B,y) as a function of wind speed is given by 


SB,y)do =pe,(B,¥)2,(B51))—-—do 3.30 
cos"(B) 
where 
z,=~tan(B)sin(y) 3.31 
z,=—tan(B)cos(y) 3.32 


and denote the surface slope where y points upwind and p(z,y) gives the fraction of a 
horizontal unit area of surface whose slopes are within z, + 1/2 5z, and z,+1/26z, 
(Basener and McCoyd, 1967). The facet slope distribution for crosswind and upwind- 
downwind situations are shown in Figures 3.4 and 3.5, respectively. The horizontal 
axis shows the value of beta, the surface normal angle from the vertical direction 
(Sidran, 1981). 

The above formulation ignores affects due to shadowing. Correction factor 
equations were developed by Wagner and also by Saunders (Sidran, 1981 and Gregoris, 
Yu, Cooper, and Milne, 1992). Omitting spectral dependencies and including Saunders’ 


correction factor S(@), the average apparent radiance from the sea is 


s(8,)t(8,) 3.33 


N,,(8768,) = 2 cos(8,) 


nr 


nr 2 
Lf [UCeFO¥, 0rd Neca -€epOpY Or ONei) 


-x9 


+(1-€F(B pb pO DN goiA8n9,)]- 


tan(6,)cos(6,,7,:9,.6/18,,7,) 


cos*(8,) }d8,dy,] 
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Figure 3.4 Facet slope distribution - crosswind case (Sidran, 1981). 
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Figure 3.5 Facet slope distribution - upwind/downwind case (Sidran, 1981). 
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Then the effective polarization is represented by 


P= * PF 3.34 


where N,,4, is the unpolarized path radiance which can be obtained with LOWTRAN 
atmospheric modeling. Fig 3.6 shows the theoretical values of polarization for 
MWIR(4-5.5um) and LWIR(8-11.5m) in and near sun glint regions. Note the negative 
values of polarization outside of the sun glint regions. Here the sea emission is greater 
than the reflected sky contribution. The negative values indicate parallel or vertical 
polarization. The effect is greater in the LWIR since the solar radiance is 
predominantly in the MWIR. (Gregoris, Yu, Cooper, Milne, 1992) 

Also, if Q represents the percent polarization of sea surface emitted radiation, 


then 


3.35 





Figure 3.7 shows emissivities and percent polarization as a function of viewing angle 8. 
Polarization is independent of surface roughness for theta less than approximately 25 
degrees. Beyond this value of theta, the polarization increases with viewing angle, but 
is partially reduced with increasing roughness. In this case an increase in polarization 
again indicates parallel polarization. Figure 3.8 shows observed data by Hale and 
Querry for emissivities and percent polarization for both reflected and emitted 
radiation for all wavelengths below 10 microns. The observed and theoretical results 


are comparable.(Sidran, 1981) 
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Figure 3.7 Emissivities and percent polarization as a function of viewing angle 
(Sidran, 1991). 
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Figure 3.8 Emissivities and percent polarization for reflected and emitted 
radiation for all wavelengths below 10 microns (Sidran, 1991). 











2. Infrared Polarization of Ship Signatures and Background Contrast 
As stated before, the infrared radiance from the sea surface is composed of two 


components. These are the thermally emitted and reflected radiation given by 


Nyp=(Ngty)!+(N goa)! 3.36 


N= P(8)Noty + 2 (8)Ngeg 3.37 


where each component depends on the reflectivity and emissivity of the sea surface 
averaged over the distribution of wave facet orientations. The emissivities for parallel 


and perpendicular polarizations are 


€1=€,c0S"(a) +e,sin*(a) 3.38 


ecg €,Sin’ (1) +€,cos”(ax). 3.39 
In terms of average apparent radiance the degree of polarization is given by 


_ <N,>-<N,> 


=v 3.40 
<N,> +<N,> 


POL 


and the radiance difference (absolute contrast) for any polarization can be written as 


AN=<Nyayget? <Nokg : 3.41 


Thus, the polarization improvement factor for contrast for horizontal filters over 


vertical filters is 


_AN,-AN, 
AN,+AN, 


3.42 


(Cooper, Lentz, Walker, Chan, 1994). 








Contrast relative to the background is most commonly defined as the difference 


between target radiance and background radiance divided by the background radiance, 


3.43 


Contrast improvement then is defined as the ratio between polarized contrast and 


unpolarized contrast, 


C, t 
Cie Ge 3.44 


Contrast between an unpolarized target against a sea background outside of 
sun glint regions should be improved with the use of polarization filters. Man-made 
targets do not indicate any predominant polarization (Cooper, Lentz, Walker, Chan, 
1994). Therefore, horizontal polarization filters would reduce the predominantly 
vertically polarized sea background radiance. Although each filter will reduce the 
unpolarized target radiance by a factor of approximately two, the reduction in sea 
background radiance should be greater. This will increase the numerator of Equation 
3.43 for horizontal polarization while reducing the denominator. Thus, an overall 


improvement in contrast would be realized. 
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IV. DATA ANALYSIS 


A. MAPTIP EXPERIMENT 


The polarization image data utilized in this thesis was obtained during the 
Marine Aerosol Properties and Thermal Imager Performance Trial (MAPTIP) held 
October 11 through November 5, 1993. The purpose of the experiment was to obtain 
atmospheric propagation data below 10 meters elevation in the 3-5m and 8-12ym 
bands. Accurate atmospheric prediction codes are necessary for properly analyzing the 
performance of electro-optical detection and tracking systems. Currently the Navy 
Atmospheric Model (NAM) used in LOWTRAN7 effectively predicts infrared 
propagation in the marine boundary layer above 10 meters. The Navy Oceanic 
Vertical Atmospheric Model (NOVAM) predicts dynamic marine aerosol composition for 
specific marine environments within the marine boundary layer. MAPTIP is expected 
to provide coastal atmospheric propagation information to extend NOVAM to the 
littoral environment (Jensen and de Leeuw, 1993). Image data used in this thesis was 
obtained during the experiment and is of the Hr. Ms. Tydemann (TDM), an 
oceanographic ship, and of a Lynx helicopter. 


B. DATA COLLECTION EQUIPMENT 


Image data was acquired during MAPTIP using the NPS AGA 780 
Thermovision Thermal Imaging system. Images were stored using AGEMA CATS data 
acquisition program. 

1. AGA 780 Thermovision Thermal Imaging System 

The version of the AGA 780 Thermovision Thermal Imaging System used 
during MAPTIP is a dual band serial scanning infrared imager. The AGA 780 
operates two scanning systems side-by-side. One system covers the 3-5.6 pm and the 
other covers 8-14 ym wavelength bands. Each utilizes a 7°x7° field of view lens that 
focuses incoming infrared radiation through a dual rotating prism system onto their 
respective detectors. The short wave lens is silicon, and its detector is an indium 
antimonide (InSb) photovoltaic detector. The long wave lens is germanium and its 
detector is a mercury-cadmium-telluride (HgCdTe) photovoltaic detector. The detectors 
are cooled to 77K by dewars filled with liquid nitrogen. The radiation passes through 


a vertically mounted prism rotating at 180 rpm followed by a horizontally mounted 
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prism rotating at 18000 rpm. The rotating prisms are synchronized electronically to 
the monitor's vertical and horizontal scans. This synchronization of the prisms serves 
to generate four fields of 100 horizontal scan lines each. Only 70 lines in each field are 
used for producing an image; the four fields are interlaced to produce an image of 280 
lines. The scanner operates at 25 fields per second. Thus, 6.25 images per second are 
generated. The detectors produce a current signal called "thermal value" proportional 
to the intensity of radiation received. This signal is amplified and converted into a 
video signal for monitor display.(AGA Operating Manual, 1980) 

2. Thermal Imager Calibration Technique 

The Stefan-Boltzmann law predicts a bolometric intensity proportional to the 
fourth power of temperature. Although the intensity of radiation received at the 
detector is presumed to be proportional to the system's "thermal values," this signal is 
a non-linear function of temperature. The AGA 780 requires calibration with a 
laboratory blackbody source. The system output of "thermal value" is recorded against 
the blackbody temperature through a range of temperatures. The calibration curve is 


generated through a least squares fit method and is of the form 


l= 





-OC 41 


where A, B, C are constants specific to the combination of scanner, lens, and filters 
used during the calibration.(AGA Operating Manual, 1980) The OC is for an offset 
correction. This correction was applied during MAPTIP to adjust for detector response 
drift (Chan, 1993). 

3. Polarizers 

The polarizers utilized during MAPTIP were Cambridge Physical Sciences 
Model IGP-228-71 with a diameter of 71 mm. The polarizing material is a 0.12 pm 
aluminum grid of 0.25 pm pitch on a substrate of barium fluoride. They have 85% 
transmissivity and 98% polarization at 3.9 pm. The polarizers were placed externally 
on the scanner lens systems and rotated manually to achieve vertical and horizontal 
polarizations.(Chan, 1993) The filters placed externally reduce the system sensitivity 
due to the added attenuation effects of the filters. Any change in system setup, such 
as including filters requires recalibration on the AGA 780 (AGA OP MANUAL, 1980). 
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4, AGEMA CATS Program 

The AGEMA CATS software allows for processing images immediately or for 
storing images to be analyzed in the future. The program is used by two PC's, each 
connected to one of the scanner channels and obtains every other field produced by the 
AGA 780. Thus CATS generates a 140x140 array for each image recorded. 


C. IMAGE PROCESSING AND ANALYSIS 


Although processing capability is available with AGEMA CATS, this thesis 
utilizes the Interactive Data Language for all image processing and analysis. 

1. Interactive Data Language (IDL) 

The Interactive Data Language (IDL) program was used to conduct the analysis 
in this thesis. This software is an array-based system with extensive capabilities in 
image processing. Numerous analysis methods and graphical representations are 
available for both programming and interactive manipulation. All images and most 
plots displayed in this thesis were produced with IDL and converted to TIFF files. The 
files were then printed with Pagemaker4. Three IDL programs were produced in the 
course of this thesis and are presented in Appendix A. Also, numerous routines were 
generated to help with the analysis of images and subsequently to assist in the 
formulation of the three analysis techniques contained in the programs of Appendix A. 

2. Image Format 

As stated earlier, all images were recorded by AGEMA CATS for processing by 
IDL. Each image is stored as a 140x140 byte file. The first 846 bytes of the file 
contain various image information tags, such as filename, date, time, offset correction 


applied, etc. An example of this "header" information is shown below. 


Thu Nov 17 17:42:44 1994 Date of ANALYSIS 

e:\ holland \oct19\lw 

filename level range offset over under date time 
OCT19100.IMG 48.80 2.00 -32.00 19600 0 10/19/93 138: 4:54 
OCT19101.IMG 44.20 2.00 -32.00 1226 13379 10/19/93 13:13:11 
OCT19102.IMG 44.20 2.00 -32.00 1175 16313 10/19/93 13:13:19 
OCT19103.IMG 44.40 5.00 -32.00 226 5716 10/19/93 13:13:29 
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The image file is read into a 140x140 array in IDL. The array values are the scanner 





generated thermal values scaled from 0 to 255. Zero represents the coldest value in 
the image and 255 represents the hottest value viewed (AGA Manual, 1980). 

3. Narcissus Effect 

All the MAPTIP images show a dark spot in the center of the image similar to 
that shown in Figure 4.1. This is known as a Narcissus spot. The Narcissus spot is 
created by retro-reflection of radiation from the detector plane by the external filters. 
Radiation from the cold detector mounting at approximately 77K reflects back toward 
the detector. The geometry of the internal optics are such that the cold detector 
radiation focuses on the detector as a small circular region near the center of the 
image. The effect of the Narcissus spot in addition to the scene obstruction 
encountered is to reduce the effective dynamic range. The extremely low values of 
irradiance from the Narcissus relative to the average scene value compresses the range 
of scene values in the upper end of the dynamic range, resulting in a loss of scene 
resolution.(Lloyd, 1975) 

The Narcissus in the system used occurs because the filters are mounted 
normal to the radiation from the detector plane collimated by the lens system. The 
reflected radiation is then focused back on the detector plane. Without the filters 
installed, the Narcissus is not present. 

Images were taken with regard to the Narcissus spot. Target images were 
placed away from the central Narcissus spot during data collection to reduce its affect. 
Although the Narcissus could in principle be removed by image processing (Chan, 
1993), it would be difficult if not impossible to replace the lost information accurately 
with other than an average or a randomly modulated average background radiance 
value. In many images, a substantial depression in background is realized that is 
difficult to ascertain from a strict viewing of the image alone. This affect can be seen 
in Figure 4.11 as not only the hole from the Narcissus, but also the depressed trough 
spanning the length of the image. 

4. Ship Location Techniques 

Calculating ship-to-background contrast requires a distinction between ship 


data and background data. The ship or other target of interest must be located within 


the image array in order to process the target data separately from the background 








data. This requires a contrast between the ship and background as given by Equation 
3.41. Three primary methods for processing the image array data are presented below. 
All three utilize a thresholding method for locating the ship. Essentially, thresholding 
means acceptance of only pixels values of the array representing brightness that are 
above a defined threshold and setting to zero all values below this threshold. Ideally, 
all ship pixel values are required to be above the threshold and all background pixels 
values below. Ship pixel values below the background level may be lost. Technically, 
to remove all the background, the threshold must be set to the maximum background 
value. Loss of ship radiance values below this value is possible. 

Care must be taken to compare appropriate background and ship values. 
Figures 4.1a and 4.2a show typical high signal-to-noise ratio images for vertical and 
horizontal polarizations, respectively. Figures 4.1b and 4.2b depict corresponding 
radiance profiles. The dark curves represent the average ship radiance profile. This is 
calculated by summing the ship radiance along each row of the image array and 
dividing by the number of ship pixels in that row. The gray curves depict the 
background radiance profile. This is calculated below the horizon by averaging all 
pixels that are not ship pixels. Above the horizon the average is calculated using only 
pixels along the edge of the image to avoid the Narcissus. Figures 4.3a and 4.4a show 
similar images for low signal-to-noise ratio images. Representative background profiles 
of the image in Figure 4.1a are depicted in Figure 4.5a and 4.6a. The upper curve in 
Figure 4.5a shows the sea radiance profile near the bottom of the image. The three 
lower curves show the sky radiance profiles. The topmost of the three curves 
represents the lowest elevation. The lower curves represent sequentially increasing 
elevation. Thus, the sky radiance decreases with elevation. Figure 4.6a shows a 
vertical slice through the image background. Analysis of this figure and Figure 4.5a 
reveals that the background varies significantly in the vertical and less in the 
horizontal ignoring effects from the Narcissus spot. The background is approximately 
constant in value below the horizon, but rises sharply above the horizon and then 
tapers off at higher elevations. Therefore, to accurately define the contrast between 
the ship and background, the ship values below the horizon should be compared with 
background values below the horizon. Likewise, ship values above the horizon should 


be compared with background values above the horizon. In addition, the pixels used to 











calculate the background average values should be as close to the ship as possible to 
represent accurately the contrast between the ship and the background in the 
immediate vicinity. 

Thresholding works well for targets with a high signal-to-noise ratio. However, 
choosing an appropriate threshold value is not clearly defined and somewhat 
subjective. Figure 4.7a shows the image from Figure 4.1a and 4.2a above using a 
threshold value equal to the background average value times 1.2. As the threshold 
value is increased to (1.3) x (background average) in Figure 4.7b and (1.5) x 
(background value) in Figure 4.7c more of the ship is deleted along the edges. From 
Figure 4.7a it is clear that choosing a threshold value is not difficult for high signal-to- 
noise ratio targets. Figure 4.8 shows the low signal-to-noise ratio image of Figures 
4.3a and 4.4a with background multiplying factors of 1.01, 1.03, and 1.05, respectively. 
In this series of images some background exists in the remaining image along with the 
ship. In this particular case, the background can be separated from calculation of ship 
radiance by confining the calculations to a region near the ship. However, in many | 
images, some background remains too close to the ship to avoid and also introduces a 
varying source of error. As the background factor is increased sequentially, thus 
increasing the threshold values, more of the background from the previous image is 
removed. However, increasing amounts of the ship image are removed, also. 
Appropriate thresholding is a balance between removing a sufficient quantity of 
background while retaining a representative amount of the target ship. This balance 
becomes increasingly difficult with lower signal-to-noise ratio targets. Often, a large 
part of the ship is not visible over the surrounding background and cannot be located 


in the image. These concepts are further described and applied in the following image 


processing programs. 
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Oct 19, 1993, (b) Average ship and average background apparent radiance profiles. Ship radiance is 
denoted by the black curve and background radiance by the gray curve. 


Figure 4.1 Vertical Polarization for high signal-to-noise 
ch 








Figure 4.2 Horizontal Polarization for high signal-to-noise ratio image (a) IDL generated image taken at 
4313 Oct 19, 1993, (b> Average ship and average background apparent radiance profiles. 
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Figure 4.3 Vertical Polarization for low signal-to-noise ratio image (a) IDL generated image taken at 0015 
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Oct 28, 1993, (6) Average ship and average background apparent radiance profiles. Ship radiance is 


= 


denoted by the black curve and background by the gray 
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Figure 4.4 Horizontal Polarization for low signal-to-noise ratio image (a) IDL generated image taken 0915 
Oct 28, 1993. (b) Average ship and average background apparent radiance profiles. Ship radiance is 
denoted by the black curve and background radiance by the gray curve. 











Figure 4.5 Horizontal profiles for high signal-to- 
noise ratio image (a) Upper curve shows sea 
background, lower three curves show sky back- 
ground, with highest of the three representing the 
lowest elevation(b) Profile through ship. Wings 
represent sea background. 





Figure 4.6 Vertical profiles for high signal-to- 
noise ratio image with horizon at position 100 (a) 
Background, (b) Profile through ship. Note cold 
spot from Narcissus where radiance drops to zero. 
Ship located at large spike centered on position 
100. 
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Figure 4.7 Thresholding high S/N ratio image. 


Left image is vertical polarization and right image 
is horizontal polarization (a) Background factor 
1.2 (b) Background factor = 1.3 (c) Bac! 
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Figure 4.8 Thresholding low S/N ratio image. 

Left image is vertical polarization and right image 
is horizontal polarization (a) Background factor = 
1.01 (b) Background factor = 1.03 {c) background 


facto 


as 


1.05. 











5. Image Processing Programs 

This section presents three processing programs developed to calculate the 
contrast improvement factor associated with the visual improvement in contrast 
apparent in the MAPTIP images. Although contrast improvement in horizontal 
polarization over vertical polarization is visually apparent in the image pairs, the 
pursuit of these programs was to quantify the contrast improvement of horizontal 
and/or vertical polarization over unpolarized images. These methods are limited to the 
apparent radiance received and recorded as image files. The calculations do not 
consider visual modulation transfer functions associated with the detection or 
identification of an image as a ship, helicopter, airplane, etc. The computations 
account for the relative differences in apparent radiance between target images and 
background. 

Crucial to these calculations is the concept of generating an unpolarized image 
from two images of opposite polarization. Stokes showed that light intensities of 
perpendicular polarizations can be added to obtain an unpolarized image or 
that the intensities of a polarized and unpolarized image can be added to obtain a 
partially polarized image (Collett, 1993). The AGA 780 response is proportional to the 
intensity of radiation received (AGA Operating Manual, 1980). Therefore, the addition 
of a vertically polarized image with the same image in horizontal polarization produces 
an unpolarized representation of the image. However, the CATS program stores the 
images scaled relative to 0-255 with zero being the lowest radiance received and 255 
the highest radiance received. Also, IDL displays images as byte array values on the 
same scale. Therefore, to produce the unpolarized image with IDL, the two polarized 


images are averaged as 


Pr aa AA 4.2 
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Image registration is another important element in the process. Since the 
images are not taken simultaneously, the ship has moved between images. Vertically 
this is usually not more than a couple of scan rows, but even that may be significant. 


However, depending upon ship speed and range, the movement horizontally may be 








substantial. This prevents the straight addition of the individual arrays. The ships 
within the arrays must be registered on to each other in order to average their values. 

a. Contiguous Area Method 

The contiguous area method was produced by W.J. Lentz at the Naval 
Postgraduate School. This program automatically finds the hottest spot in the image 
and works through a contiguous set of pixels to the edges of the ship where the 
temperature drops. The spot may be chosen if the hottest spot is background clutter. 
Pixel values above a set background threshold value are chosen and all others set to 
zero. The background threshold value may be varied. Each horizontal line in the 
image is ascribed a background temperature based on a smoothed average along the 
edges of the image. The values are used to find variations in temperature along the 
horizontal lines that could represent the ship. A small temperature difference might 
miss much of the ship. For this reason the search is iterative over the entire image 
using an edge threshold. Once the ship is detected, another threshold is taken in a box 
around the ship. This corrects for effects such as the Narcissus spot superimposed on 
the image. The output is a file that gives the contrast of the ship to background in a 
local box that extends to the limits of the ship.(Lentz, 1994) 

The thresholding idea was first used to discern the ship from the 
background in a joint experiment with NOSC. In that experiment, the AGA 780 was 
used without polarizing filters, and it was relatively easy to find the ship. In the 
MAPTIP experiment, the signal-to-noise ratio was much worse both because of the 
polarizers that reduced the signal but also because the targets were generally much 
further away. The earlier program either could not find the ship or would include too 
much background with the ship. In addition, the sink-hole effect in the center of the 
AGA pictures caused by the Narcissus effect produced additional distortions that 
caused problems. The essential idea to correct these problems is the contiguous point 
finding program in IDL called search2D used in the subroutine shipfind.pro. It starts 
from the designated hot spot and eliminates noncontiguous points from the picture. 
Thus what remains is a continuous matrix of pixels values that are above threshold. 
All pixels not connected to the hot spot are set to zero. Then a box is defined at the x 
and y limits of the ship. Then the threshold which had been at the edge of the picture 
is redefined for the vertical edges of the box in which the ship is contained. The 
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process is repeated producing a cleaner ship with a lower allowable temperature 
difference for background.(Lentz, 1994) 

Since ship pictures are taken in pairs, the program takes the next image 
and compares it to the previous image and combines the ship images to produce an 
unpolarized ship image. The ship may have moved slightly from picture to picture, so 
the program adjusts for the motion of the box. One polarization or the other will have 
an improved contrast over the background resulting in a larger ship image with more 
pixels. In some cases, for example, the ship will consist only of a smokestack with one 
polarization and in the other the entire ship will be clearly defined. In order to 
properly compare the contrasts, the ship with the larger number of pixels is used to 
redefine the ship outline in the poorer picture. The box and ship sizes are the same for 
both pictures giving a better calculation of the polarization improvement. The ship 
contrast is formed between all the points on the ship and all the points in the box 
containing the ship that are not the ship. This is done for both polarizations and for 
the sum of both polarizations. (Lentz, 1994) 

The significant advantage of this program is its autonomy. Many image 
pairs can be processed in a minimum amount of time. However, this program does not 
compare the entire ship as observed. Since only contiguous areas are used, only the 
contiguous area containing the hot spot of each ship image is used to calculate 
contrast. b. Outline Method 

The “outline method" is an attempt to overcome the difficulty of finding 
a low signal to noise ratio ship image and is provided in Appendix A as the program 
outlin15.pro. The high signal-to-noise ratio ship image of Figure 4.1 was used to 
generate a block outline of the beam aspect of the ship. The outline was referenced to 
an IDL generated box function. This allowed for processing ship targets at various 
ranges since at longer ranges the ship takes up less area in the field of view. When 
the box is placed around the ship touching only the stack, bow, stern, and waterline, 
the outline is scaled to fit the ship. If only part of the ship is visible, then judgment is 
required to align the box properly. 











Figure 4.9 IDL ship image with outline in place. 
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Once the box is placed around the ship, average background values are 
calculated row by row using five pixels to either side of the box. The outline is drawn 
referenced to the box as shown in Figure 4.9 All processing was performed on raw 
data image files. Figure 9 is an expanded file of dimensions 280 x 280. This 
expansion process consists of doubling the size of each side of the image array. The 
new pixel positions are filled by linearly interpolating between the original pixels. 
Note the "steps" created during the expansion along the bow and stern outlines. All 
pixels outside of the box are then set to zero. The average background value of each 
row multiplied by an input factor to account for pixels above the average is used to set 
to zero any remaining pixels within the outline below the threshold. Originally, the 
background average values were calculated using the five pixels on either side of the 
outline. Ship edge effects and even a slight error in placement of the IDL box lead to 
significant errors in background average values when part of the ship registers as one 
or more of the 10 background pixels, i.e. part of the ship falls outside of the identified 
boundary. Therefore, five pixels on either side of the box were used to generate 
average background values for each row. 

At this point all that remains in the array is the visible portion of the 
ship. Then the program steps through each row calculating the average ship value of 
that row using the remaining pixels in that row. This produces a profile of the average 
ship values which can be displayed graphically as in Figure 4.1b and 4.2b. A second 
image and profiles are produced for the perpendicular polarization. The profiles of 
contrast and contrast improvement may also be represented similarly. Then using 
these values average contrast and average contrast improvement for above and below 
the horizon may be determined. 

However, coupling of the ship areas is required. Calculating average 
ship values of the individual polarizations without coupling the areas is comparable to 
comparing images of two different targets. However, this method erroneously 
accounted for this coupling. First, the "outline method” assumes that the ship images 
are registered exactly due to the outline placement. Therefore, each row of one image 
is compared to its respective row in the orthogonal polarization. A 140 x 140 image 
window on the computer monitor is 38mm x 38mm. This requires accurate placement 


of the box around the ship to within 0.27mm. The assumption of proper registration is 
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invalid. Second, coupling was achieved by comparing the number of ship pixels in a 
given row to the number of ship pixels in the same row of orthogonal polarization. The 
row with the most pixels and thus the larger area must represent the minimum actual 
ship area along that elevation. Then the sum of pixels along both rows was divided by 
the larger number of pixels to obtain the effective average ship radiance at that row. 
This erroneously reduces the apparent average radiance of the row with the lower 
number of pixels. 

Several problems exist with the "outline method" in addition to the 
coupling described above. It is difficult to place the box properly around low signal-to- 
noise ratio ship images. Improper placement and sizing results in the outline pattern 
not including all of the ship on one side and including excess background on the 
opposite side. The major drawback with the "outline method" is that it is restricted to 
beam aspect only. The image array may be inverted without loss of information in 
order to utilize the outline on either the port or starboard aspect. However, more 
detailed programming is required to allow for other than beam aspects although this 
method is probably too cumbersome to provide accurate outlining in other than beam 
aspects. 

c. Box Method 

The "box method" is the most simple of the three methods presented. In 
fact it was the first one attempted, but later dropped to permit experimentation with 
the previous methods. This methcd simply places a box around the ship and 
determines the average background value about the ship row by row using the five 
pixels to either side of the box. All pixels outside of the box are set to zero. Then a 
threshold is obtained using the average background value and an additional amount to 
account for background values above the average. As before the values below this 
threshold within the box are set to zero. However, this method contains a significant 
change in the fundamental interpretation of the calculation of average ship value. 

The "outline method" calculated the number of pixels per row that were 
ship pixels for both vertical and horizontal polarizations. The larger number of ship 
pixels along with the sum of ship values per row were used to calculate average ship 
value for vertical and horizontal polarizations. The two polarizations were then 


averaged to produce an unpolarized image. This method was flawed because it 
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improperly couples together the two images with respect to ship area. The image with 
the smaller number of pixels is considered to have an area greater than the actual 
area. Thus pixel values of zero are included in calculating its average radiance. Since 
the ship pixel values must be at least that of the background, assigning them a value 
of zero is incorrect. This error reduces the calculated average radiance. The "box 
method" corrects for this flaw. 

6. Box Method Generation of Unpolarized Image 

In the "box method" the respective ship areas in the vertical and horizontal 
polarizations are coupled. The ship areas are represented by the number of ship pixels 
in the image. Ship area is subdivided further with respect to ship area above and 
below the horizon. The image with the larger area or greater number of pixels above 
or below the horizon defines the ship area for that region. In the other image of 
smaller area, the "missing" pixels are given a value equal to the background value. 
The sum of all pixels in their respective regions is divided by the total number of ship 
pixels in that region. Now both ship images contain the same area of ship surface. 
This couples the ship areas by referencing the average ship values to the same ship 
area. The motivation for this is that while part of the ship in one image cannot be 
seen over the background it actually is being seen with the value of the background. 
Now the average ship value of each polarization is being referenced to the same size or 
area of ship although this area is not necessarily the true ship projected area. Even 
the image with the larger area may not represent the total projected area if part of its 
radiance is the same as the background average. Without this coupling provided, the 
smaller ship image average value is calculated erroneously lower than the actual value 
leading to erroneously low contrast and contrast improvement values. In certain 
instances, the "outline method" produced negative contrast values which was 
impossible since the detectable ship is of higher radiance than the background. The 


negative contrast values indicated the existence of and led to the identification of the 


area coupling flaw. 











Some image information is lost with the "box method." Since the average ship 
values above and below the horizon are calculated as the sum of all ship radiance 
values divided by the total number of pixels, the row by row profile information is lost. 
If the two ship images could be accurately registered, then a row by row comparison 
would be possible. However, since the relative radiance values of the ship images with 
respect to polarization do not differ beyond what could be attributed to noise, the 
profile information from these images would provide no new meaningful information. 

7. Ship Signature Polarization 

In an attempt to analyze the polarization differences between horizontal and 
vertical polarizations, two methods were used. The first is a three dimensional surface 
plot or surface representation. The second is a cross-sectional view of the image array 
or profile representation. 

a. Surface Representation 

A three dimensional surface plot of an image array is shown in Figure 
4, 10. The top image is of the vertical polarization of a ship image and the bottom is of 
horizontal polarization. Figure 4.11 shows the same plot rotated about the vertical 
axis. The effect of the Narcissus can be seen better from this angle. The Narcissus 
can be seen as the depression or hole in the middle of the plot. Additional damage by 
Narcissus depression can be seen emanating away from the primary "hole" in the form 
of a trough. This cold spot reduces the effective dynamic range of the array values by 
forcing the target radiance values to be compressed at the upper end of the dynamic 
range. 

Comparison of the two surface plots shows that it is difficult to discern 
ship radiance value differences. What differences are present cannot be justified as 
anything other than noise. Also, it is difficult to determine polarization affects on the 
background. Lastly, this data provides no numerical information, although it does 
provide some visual insight. Some images such as the Lynx helicopter in Figure 4.12 
exhibit flat surfaces in these plots as shown in Figure 4.13. These flat surfaces 
represent pixels values that exceed 255. Their actual apparent values are 


undetermined. 
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b. Profile Representation 

An alternate representation is shown in Figure 4.5b and 4.6b. These 
show horizontal and vertical slices through the ship. All the slices of either vertical or 
horizontal profiles can be averaged to obtain a representative profile of the ship. The 
Narcissus spot and the nature of the background lead to averaging the background 
pixel value of each row instead of each column to produce a representative vertical 
profile. Vertical profiles complement the process of comparing ship row values with 
their respective background row values. Comparison of vertical profiles of horizontal 
and vertical polarizations again shows little polarization of radiation from the ship. 
However, vertical profiles show the effect polarization filters have on the below-horizon 
background radiation received. This is demonstrated in Figures 4.14 and 4.15. These 
graphs show superpositions of the average horizontal and vertical background radiance 
values and ship values. The unpolarized radiance would be the average of the two and 
positioned midway between the lines shown. The background radiation received below 
the horizon is significantly reduced with the horizontal polarization filter. A small 
difference in radiation received exists above the horizon. Thus, the sea background 
radiation is predominantly vertical polarized. No discernable difference exists along 
the ships in Figures 4.14b and 4.15b and thus the ship radiance is essentially 
unpolarized. 

e. Average Value Calculation 

The average value of ship radiance for a given polarization may be 
calculated with two methods. The first method calculates the ship average values for 
both above and below the horizon. The "box method" program reduces the current 
image array to only ship pixels and all other pixels are set to zero. The horizon is 
marked with the mouse cursor. The average ship value above the horizon is calculated 
simply by summing the pixels above the marked horizon and dividing by the number of 
ship pixels above the horizon. The average ship value below the horizon is calculated 


similarly. 











Figure 4.10 Surface plot of image in Figure 4.1a. Flatter base region repre- 
sents background radiance value. Elevated portions show ship radiance 
values (a) Vertical polarization (b) Horizontal polarization. 











Figure 4.11 Surface plot of Figure 4.10 rotated (a) Vertical polarization (b) 
Horizontal polarization. 
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Figure 4.12 IDL generated image of Lynx helicopter 
(a) Vertical polarization (>) Horizontal polarization. 
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Figure 4.13 Surface plot of helicopter in Figure 4.12. Large spike shows the engine 
heating. Smaller peaks are from local hot spots on the helicopter (a) Vertical polarization 
(b) Horizontal polarization. 














Figure 4.15 Characteristic apparent radiance 
profile comparisons for low S/N ratio image (a) 
Background comparison (b) Ship comparison. 


Figure 4.14 Characteristic apparent radiance 
profile comparisons for high S/N ratio image (a) 
Background comparison (b) Ship comparison. 
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Alternatively, the average value of ship radiance for a given polarization 
is calculated for each row. This average is calculated simply by summing each row 
between the top and bottom of the box and dividing by the number of ship pixels in a 
row. These averages then constitute the image representative vertical ship profile. 

8. Background Polarization 

Background polarization effects were analyzed using image profiles. Significant 
vertical polarization was observed below the horizon. Less vertical polarization was 
observed above the horizon. 

a. Average Value Calculation 

The average value of background radiation for a given polarization is 
calculated for each row. The procedure varies with the location of the target with 
respect to the Narcissus spot. For images where the target is above the Narcissus 
spot, the background above the ship (i.e. above the box) is calculated simply by 
summing each row and dividing by the number of pixels in a row (i.e. 140). Below the 
ship the first and last 20 pixels in each row are summed and divided by 40 to 
determine the average background value below the ship. Along the ship the five pixels 
to either side of the box are averaged to find the average background value along the 
ship. These averages then constitute the image representative vertical background 
profile. 

b. Characteristic Profile 

Previously, Figure 4.5 and 4.6 showed both horizontal and vertical 
profiles of an image. The horizontal profiles depict a background that is essentially 
constant at a given elevation disregarding the affect of the Narcissus. The vertical 
profiles depict a widely varying background with respect to elevation. A target is 
discernable from the background in the immediate vicinity if there exists an 
appreciable contrast between the ship radiance and the background radiance. The key 
is "immediate vicinity." Contrast must be based upon appropriately related ship and 
background values. It is not appropriate to compare ship radiance values below the 
horizon with background above the horizon. Therefore, the ship radiance values along 
a given row must be compared with background radiance values along that same row. 


Thus, all the pixels in each row are averaged to generate a characteristic profile. 











The characteristic profiles of Figures 4.14a and 4.15a show a significant 
reduction in background radiance below the horizon while no significant difference 
exists in the radiation from the target ship. It is this reduction in background 
radiation that polarization filters can take advantage of and increase target to 
background contrast. By aligning the filters to block the predominantly vertically 
polarized background radiation the contrast will be enhanced. 

9. Ship to Background Contrast 

Once the ship and background average values have been determined, the 
contrast values may be calculated. Contrast is defined as a difference in radiance 
value between the target and surrounding background. This difference may be positive 
or negative depending on whether the target is of greater or lesser radiance than the 
background. Equation 4.3 below provides the mathematical definition of contrast. 
Ideally, contrast results from any one or a few adjacent pixels being of greater or lesser 
value than the surrounding background pixels. This introduces the concept of clutt=r. 
Clouds may reflect radiation into the detector with a value comparable to the ship 
radiance. The ship targets contain more than one individual pixel as is necessary to 
achieve classification or identification. So there exists a question as to what is 
detection that is not addressed by this thesis. As stated before, this thesis does not 
include modulation transfer functions for determining detection. 

10. Contrast Improvement 

The quantities of interest are the improvement in vertical or horizontal 
polarization contrast over the unpolarized contrast. Provided there is a significant 
improvement in one of the polarizations, then such a polarizer may prove useful in an 
actual target detection or tracking system. It is reasonable that a variable plane of 
polarization filter might be included in system design. 

a. Definition of Contrast Improvement 


Target to background contrast can be defined as the normalized 


radiance difference between target and background where N,,, and N,,4 are apparent 











radiance values, respectively. Therefore, the contrast for vertical, horizontal, and 


unpolarized target images is 
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Then the contrast improvement for vertical and horizontal polarizations is defined as 
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b. Background Dependency 

As stated before, contrast calculations are dependent upon the 
background utilized in the calculations. If the background average is determined using 
the entire range of values from the lowest radiance at high elevations to the highest 
radiance just above the horizon, then this average may be significantly higher or lower 
than the average background value in the vicinity of the target. Care must be taken to 
associate the appropriate background with the target. This is most effectively done by 
comparing the below horizon background with the ship component below the horizon 
and the above horizon background with the ship component above the horizon. A more 
accurate method would be to do a row by row comparison. This is easily done with one 
image of vertical or horizontal polarization. However, for meaningful analysis the 


images must be summed. This leads to the problem of registration covered earlier. 





c. Ship Area Reference Dependency 

The contrast calculations depend on the area coupling discussed 
previously. Originally, row by row analysis was performed. This entailed summing 
the ship pixels along each row. The same size box was used in each polarization and 
assumed to register the ship image properly. The same row of each image was 
compared by using the larger number of ship pixels in each row from the respective 
vertical and horizontal images. This often leads to negative contrast values, since the 
calculated N,,, may be erroneously low. This problem was resolved by first 
acknowledging the error in ship registration. Then instead of a row by row analysis 
the ship radiance was calculated as an average either above or below the horizon. 
Also, the difference in ship pixels was compensated for by assigning the background 


average to the "missing" pixels. 


D. RESULTS OF ANALYSIS 


Contrast improvement calculations were generated by the "contiguous area” 
and "box" methods. Visual review of the MAPTIP image data reveals a greater ship to 
background contrast in the horizontally polarized images than in the vertically 
polarized images. Of interest is the degree of contrast improvement achieved with 
polarizers over that without polarization filters. Unpolarized image data is not 
available for these images and must be artificially constructed. Both the contiguous 
area and box methods generate unpolarized images and calculate the contrast 
improvement obtained with vertical and horizontal filters as described above. The 
results of these calculations are provided in Appendix B. 

Both methods show con:’ ast improvement in the horizontal polarization. 
Sometimes even vertical polarization filtering shows a contrast improvement over 
unpolarized -nages. Horizontal polarization consistently provides a greater contrast 
improvement over vertical polarization except within a strong sunglint channel. 

Tables 4.1 through 4.5 show results from selected images throughout different 
times of the day. Table 4.1 images were recorded in the early morning, Table 4.2 in 
mid-morning, Table 4.3 in mid-day, Table 4.4 in mid-afternoon, and Table 4.5 in the 
evening. Each table is in two parts. Part (a) lists the average radiance values for 
target and background both above and below the horizon as calculated from the "box 


method." Part (b) lists contrast and contrast improvement values. The last three 
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columns depict polarized target-to-background contrast (C), unpolarized target-to- 
background contrast (C,,), and contrast improvement (C,,,,) as calculated from the 
“contiguous area method." 

These values of contrast improvement require qualification. The results of both 
methods are not strictly comparable. Errors are inherent in both methods due to 
misidentification of the horizon or ignoring horizon effects. The contiguous area 
method calculates a ship contrast irrespective of the horizon. This is not a significant 
miscalculation in most low signal-to-noise ratio images. In these cases, the horizon is 
seldom visible and the apparent radiance just above and below the horizon is similar. 
This can be seen from the example low signal-to-noise ratio image of Figure 4.4. The 
results of calculations for this image are tabulated in Table 4.2 as file OC28L025 and 
OC28L026. The "box method" contrast and contrast improvement below the horizon is 
nearly the same as the overall contrast and contrast improvement calculated from the 
"contiguous area method." However, in images where the horizon difference is large, 
this error increases substantially. Image files OC19L105 and OC19L106 of Table 4.3 
show the results for the example high signal-to-noise ratio image of Figure 4.1. Since 
in such images the sky radiance just above the horizon is much larger than below the 
horizon, then the overall ship to background contrast is falsely reduced because the 
average background is erroneously increased by averaging above and below 
background values. This results in the negative contrast values of several image pairs 
in the "contiguous area" method data as seen in Table 4.2b for file OC25L084. The 
"box method" avoids this error. However, accurate identification of the horizon is 
difficult. This also leads to an error in calculation of the above and below horizon 
background values. 

Table 4.6 compares the "box method" contrast improvement (C,) for horizontal 
polarization below the horizon and "contiguous area method" for contrast improvement 
(C,) for horizontal polarization to two previous efforts to calculate contrast. The 
contrast in C, was calculated without correctly coupling the ship areas. In this method 
the ship image of smaller area was assumed to have the same area as the larger ship 
image. The total pixel radiance value was divided by this larger area which means the 
missing pixels were set to zero. This reduces the average ship radiance and contrast 


for that polarization. This reduced ship radiance also reduces the unpolarized ship 
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radiance and contrast which is calculated from both polarized radiance values. 
Therefore, opposite polarization contrast improvement will be increased when its 
contrast is divided by the reduced unpolarized contrast. If the vertically polarized 
image is of less area, then the horizontal contrast improvement will be inflated as 
shown in Table 4.6 as compared to C, which sets the missing pixels values of the 
smaller image equal to the background radiance. The last column shows ship-to- 
background contrast calculated from mean radiance data shown by Chan (Chan, 1993). 
This mean radiance data consists of the average ship or background radiance within a 
box placed either within the ship target or a box containing background only. 

The "contiguous area" method is semi-autonomous. Once a threshold value is 
input, the only operator action necessary occurs when the program experiences a 
conflict in location of the ship hot spot. This spot is identified by the operator and the 
program continues. In several instances, the "contiguous area" method includes a 
large and significant portion of background in its calculation of ship average radiance 
value. The "box method" also avoids this by allowing the operator to set the size of the 
box to a tight fit around the ship. 

Also, the "contiguous area" method occasionally does not include all of the ship 
in calculating average radiance. This happens when portions of the ship in both 
horizontal and vertical polarizations are in separate blocks of pixels. For, example 
some images will have the ship stack appear separated from a larger area of ship hull. 
The program will only use one of these areas in the calculation. For strict contrast, 
the stack area which is of a greater apparent temperature should be chosen. This 
leads again to the concept of detection. If the stack is to be chosen only because of 
higher apparent temperature, then in all images only the highest ship pixel value 
should be compared to the surrounding pixels. For the stack, the surrounding pixels 
are other ship pixels. The concept of pixel to pixel contrast works if the target is one 
pixel in area. The "contiguous area" method compares the ship "hot chunk" to 
background whereas the "box" method compares the entire infrared visible ship to 
background. The former applies more to pixel detection and the latter applies more to 
target detection. 

Apparent in the "box" method data is the greater contrast improvement in 


below the horizon compared to that above the horizon in horizontal polarization. This 
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results from the sea background emission which is predominantly vertically polarized 
(Cooper, Crittenden, Milne, Walker, Moss, and Gregoris, 1992). The sky background 
immediately above the horizon is predominantly unpolarized where there exists a 
significant horizon contrast as depicted in Figure 14a. This figure shows no 
measurable change in apparent radiance in either polarization. This can only occur if 
the radiance is unpolarized. 

Definitions of the results in Tables 4.1 through 4.5 are given below. The 


following values were calculated by the "box method." 


(N).1y = apparent radiance above the horizon for given polarization 

(N)-) = apparent radiance below the horizon for given polarization 
(N,,)aw = apparent radiance above the horizon for unpolarized image 
(Nypeei = apparent radiance below the horizon for unpolarized image 
C.v = contrast above the horizon for given polarization 

C,.; = contrast below the horizon for given polarization 

(C,,)aby = contrast above the horizon for unpolarized image 

(C,,)s = contrast below the horizon for unpolarized image 

(Cimp)abv = contrast improvement above the horizon for given polarization 


(Cinp)be) = contrast improvement below the horizon for given polarization 
The following values were calculated by the "contiguous area method." 


= contrast for given polarization without regard for horizon 


u 


C 
C,, = contrast for unpolarized image without regard for horizon 
C, 


= contrast improvement without regard for horizon 


ip 
imp 
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vertical/horizontal image pairs and unpolarized image witn respect to target or 
background source and horizon. Images taken during the early morning. "V" and "H" denote vertical and horizontal polarizations, 
respectively. 
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horizontal polarization as_ calculated by "box method." Last three columns show overall contrast and overall contrast improvement as 
calculated by "contiguous area method." Images taken during the early morning. "V" and "H" denote vertical and horizontal 
polarizations, respectively. 
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ple 4.2a Radi verticai/horizontal image pairs and unpolarized image wit respect to target or 
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respectively. 
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able 4.2b. The first six columns show contrast and contrast improvement values for selected image pairs with respect to horizon for 


horizontal polarization as calculated by "box method." Last three columns show overall contrast and overall contrast improvement as 


calculated by “contiguous area method." Images taken during the mid-morning. "V" and "H" denote vertical and horizontal 
polarizations, respectively. 
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able 4.3a Radiance values for seiected vertical/horizontal image pairs and unpolarized image witn respect to target or 


background source and horizon. Images taken during the mid-morning. "V" and "H" denote vertical and horizontal polarizations, 
respectively. 
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Table 4.3b. The first six columns show contrast and contrast improvement values for selected image pairs with respect to horizon for 
horizontal polarization as calculated by "box method." Last three columns show overall contrast and overall contrast improvement as 
calculated by "contiguous area method." Images taken during the mid-morning. "V" and "H" denote vertical and horizontal 


polarizations, respectively. 
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Table 4.4a Radiance values for selected vertical/horizontal image pairs and unpolarized image wit respect to target or 


background source and horizon. Images taken during the mid-morning. "V" and "H" denote vertical and horizontal polarizations, 
respectively. 
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Table 4.4b. The first six columns show contrast and contrast improvement values for selected image pairs with respect to horizon for 
horizontal polarization as_ calculated by "box method." Last three columns show overall contrast and overall contrast improvement as 
calculated by "contiguous area method." Images taken during the mid-morning. "V" and "H" denote vertical and horizontal 
polarizations, respectively. 
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contrast improvement values for selected image pairs with respect to horizon for 


horizontal polarization as calculated by "box method." Last three columns show overall contrast and overall contrast improvement as 
calculated by "contiguous area method.” Images taken during the mid-morning. "V" and "H" denote vertical and horizontal 
polarizations, respectively. 





OCI9LIV/12 
OC19L122/28 
OC19L124/25 


Table 4.6 Contrast improvement comparison. C, calculated without setting "missing" 
pixels in smaller image to background value, C, saleulated "box method", C, calculated 
using "contiguous area method," C, calculated using mean radiance data obtained by 
Chan. C, and C, represent contrast improvements below the horizon for horizontally 
polarized images. C, and C, are contrast improvements for horizontal polarizations. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The below horizon contrast improvement values confirm that a significant 
improvement in detection ability is possible with horizontal filtering. The below 
horizon values are significant because in most operational environments a sea 
background will dominate target scenes. For an aircraft flying as low as 500 feet, the 


distance to the horizon is 25 nautical miles given by 


D=1.14y(h) 5.1 


where D = distance to horizon in nautical miles 
h = height of observer in feet 
(Tobin, 1974). Then the distance to a 150 ft high surface vessel that is completely 


below the horizon is 22.5 nautical miles as shown below. 





Figure 5.1 Observer to horizon geometry. 


The surface is approximated to be flat in the diagram. The target ship height is 
denoted by TS, the observer height is denoted by OP, and the horizon by H. FLIR 
acquisition ranges are often displayed as a function of target-to-background 
temperature difference. For a typical 400 ft long destroyer from a beam aspect, 


unclassified acquisition ranges show a maximum acquisition range of less than 20 








nautical miles (Cooper class notes, 1994). Thus a significant operational improvement 
is possible with infrared detection and imaging systems using horizontal polarizing 
filters because the filters would reduce sea component of background radiance. 

These calculations are for detection only. The "box" method more closely 
approximates recognition improvement by using the entire ship radiance in the 
average radiance calculation. Further consideration needs to be applied to the 
resolution requirements necessary to recognize a detected target. Target detection in 
the presence of clutter is felt by some to require some resolution capability - e.g. FLIR 
resolution requirements are established at 1.5 line pairs for detection (Rosell and 
Harvey, 1979). Given a system resolution capability requirement, the improvements in 
detection and recognition of horizontal filtering over no filtering should be explored. 
As an example, some image pairs show a sizeable increase in total ship area detected 
with horizontal polarization (see Figure 4.9). A system's resolution requirements 
would be less with the increased area of ship radiance received in excess of background 
radiance received. 

It is difficult to develop a technique for evaluating the contrast which would be 
observed and be applicable to detection/recognition criterion. As a minimum the 
method chosen must be consistent with visual results. The "box method" is believed to 
be such a method. Validation of that conclusion is needed. Contrast improvement 
refers only to apparent contrast. This must be dependent on the range at which the 
target is viewed. Range dependence of the contrast should be included and validated. 

Image registration remains a problem. The inability to properly sum two 
images precludes correct generation of an unpolarized image. A related problem is the 
time difference between image pairs. The atmosphere is dynamic and thus the 
radiance received is time dependent. The atmospheric changes between images are 
unknown and may or may not be significant. 

One possible solution to both problems is in progress. This involves a lens 
adapter system to be mounted on the AGA 780. This system serves two purposes. 
First, the design is such that the Narcissus spot will be removed. Second, the system 
splits the incoming beam of radiation. Each beam then passes through either a 
horizontal or vertical polarizing filter. The resulting image will contain both 


polarizations each 140 x 70. Since these images were taken simultaneously, then they 
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will be registered perfectly. Then they may be summed to produce an unpolarized 
image. Once this system is constructed and mounted on the AGA 780 it must be 
calibrated and tested to verify or disprove the assertion that an unpolarized image may 
be produced with the AGA 780 by summing a vertical and horizontal image. 
(Skowroneck, 1994) 

Environmental data acquired from the MAPTIP experiment could be used for 
further analysis of the MAPTIP images. Calculations in this thesis were based on 
apparent radiance. This is the radiance actually received at the AGA 780. This is not 
the radiance emitted from the target. Contrast and contrast enhancement are relative 
- based on relative radiances. This does not pose a problem as would absolute 
radiances. With the environmental data from MAPTIP, an atmospheric transmission 
model such as modified LOWTRAN6 may be used to determine the actual radiance of 
the ship and background. This information is necessary for proper modelling of 


electro-optic system requirements. 











APPENDIX A. 


A. OUTLINE METHOD PROGRAM 


;outlin14 constructs outline of ship referenced to box placed around ship : 
box must touch top, bow tip, stern tip, waterline 


lorder=0 

print,'flirs makes movie of flirs' 

close,1 

;close in case of crash 

intfn=0 

shpnmbr=1 

ct=0 

loadct,ct 
:window,0,xsize=280,ysize=280,xpos=0,ypos=0 
on_ioerror,problem 

path1=' . 

file1=' : 

print,'file size must be 20446 with header' 
read,'input path to *.img files',path1 

if path1 ne " then path=path1 

print,'path= ',path 

read,'input type of files eg *.img or names to search for',filel 


if filel eq " then filel='*.img 
files=findfile(path+'\'+filel ,count=num) 
for i=0,num-1 do print,files(1),i 

loop: print,n_elements(files)-1,' last number" 
read,'input i for number to plot or -1 to stop’,i 
if i lt 0 then goto, problem 

istart=i 

read,'input stop number',istop 

for i=istart, istop do begin 

f=files(i) 

openr,1,f 

print,'opening ',f,i 

head=bytarr(846) 

y=bytarr(140,70) 

yl=bytarr(140,70) 

yy=bytarr(140,140) 

readu,1l,head 

readu,1,y 

readu,1l,y1 

close,1 

index=indgen(70) 
yy(*,index*2)=y1(* index) 
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yy(*,index*2+1)=y(*,index) 
if(max(yy) eq 255) then begin 
yl=yy( where(yy gt 254)) 
nover=n_elements(y1) 
endif else begin 
nover=0 
endelse 
yy=yy(*,139-indgen(140)) 


0 oR ae a a a AEE a ae a AE A EE EE OK OE EER OE KE EE 
’ 


;save ship images 
FA AOELELES ES ELLER AEL ETE PREM O REED ES TT 
shipl=bytarr(140,140) 
ship2=bytarr(140,140) 
if shpnmbr eq 1 then begin 
shipl=yy 
shipll=yy 
shpnmbr=2 
titlel='vert' 
window, 1,xpos=0,ypos=0,xsize=140,ysize=140,title=title1 
tvscl,ship1 
endif else begin 
ship2=yy 
title2="hor' 
window,2,xpos=0,ypos=160,xsize=140,ysize=140,title=title2 
tvscl,ship2 
shpnmbr=1 
endelse 


endfor 


eA EK RR ER RR EB CE RE EE EE EK 
’ 


sinvert outline for port beam aspect 
EEEREOE EMEP ELE TINE S TS TREE SED EU Te 
read,'invert outline for port beam aspect? yes=1 no=0' aspect 
if aspect eq 1 then begin 
ship11=ship11(139-indgen(140),*) 
ship2=ship2(139-indgen(140),*) 
wset,1 
tvscl,ship11 “i 
wset,2 
tvscl,ship2 
endif 


SHEE NVvert signal rere 
yl=1 


y=1 
heRRRRE HOS Ot variables******#*** 
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yy=congrid(yy,140*2,140*2, Interp=1) 
slorder=1 

loadct,0 

swset,0 

print,’Scaled plot of data in 0' 

stvscl,yy 

range=double(head,770) 
level=double(head,762) 
dlevel=double(head,618) 
emiss=double(head,626) 
ct=double(head,650) 
r=double(head,570) 
b=double(head,578) 
f=double(head,586) 
amt=double(head,674) 
at=double(head,666) 

pl=emiss*ct 

p2=(1-emiss)/emiss 

p3=1.0-ct 

iamb=r/(exp(b/amt)-f) 
iatm=r/(exp(b/at)-f) 
y=range*(float(yy)-128.)/256.0+level+dlevel 
y=y/p1-p2*iamb-p3/p1*iatm 
t=b/alog(abs((r/y+f)))-273.15 

; get rid of old variables 
s;delvar,y,yl,yy,head 


shpavel=fitarr(140) 
shpavel(*)=0.0 
shpave2=fltarr(140) 
shpave2(*)=0.0 
shptotl1=fltarr(140) 
shptot1(*)=0.0 
shptot2=fitarr(140) 
shptot2(*)=0.0 
sindex1=fltarr(140) 
sindex1(*)=0 
sindex2=fltarr(140) 
sindex2(*)=0 
shppix1=fltarr(140) 
shppix1(*)=0.0 
shppix2=fltarr(140) 
shppix2(*)=0.0 
outline=fltarr(140,140) 
outline(*,*)=0.0 
outline1=fitarr(140,140) 
outline1(*,*)=0.0 
outline2=fltarr(140,140) 








outline2(*,*)=0.0 


redo:read,'enter BF value',dd 
BF=dd 

image=1 

nextimg: 

wset,image 
zz=boxfun(x0,y0,nx,ny,init=intfn,/message) 
intfn=1 
box=tvrd(x0,y0,nx,ny) 
bkgdave=fitarr(140) 
bkgdave(*)=0.0 

delta1=0 

delta2=1 


eR A ea a EA a a a EK a A KE ACE OR A A EK AK IE RE EE RR KEKE 


> 

sbox or outline 

eH EA HCI KE KK EK A KE I KO A KK AK KK RR RE RE EK AR EE EEK KE 
> 


OEE AA a IE CO EA A OK HE EK A A OK OK A AE A A a EE AC A I A BE BR EE 
’ 


sassign image array to outline array and mark horizon 
SERED EEEANECERT ESSA UE SELES ESTE EAE ROR TERETE RE TO Te ee 
if image eq 1 then begin 
outline=ship11 
boat=ship11 
print,'click on horizon or bottom to exit’ 
replotl: cursor,xv,yv,3,/device 
if yv It 5 then goto,exit1 
xvl=xv 
yvl=yv 
print,'horizon=',yv1 
print,'click on point or bottom to exit’ 
goto,replot1 
exit1:print,'exiting' 
endif 
if image eq 2 then begin 
outline=ship2 
boat=ship2 
print,'click on horizon or bottom to exit’ 
replot2: cursor,xv,yv,3,/device 
if yv lt 5 then goto,exit2 
XV2=xV 
yv2=yv 
print, horizon="',yv2 
print,'click on point or bottom to exit’ 
goto,replot2 
exit2:print, ‘exiting’ 
endif 


eR KAA AR RK RR RE BR A ER EK A ER EB KEKE EEE 
? 
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;calculate background average values 
ERC TELR DRE RTE Ree ROC Nee NN SEE TOT eee ee ea TT Tee 
if yO+ny/2.0 It 69 then begin 
sbox below narcissus 
sbackground average below box 
for n=0,y0-1 do begin 
bkgdave(n)=total(outline(0:139,n))/140.0 
endfor 
sbackground average around box 
for n=y0,y0+ny-1 do begin 
bkgdave(n)=(total(outline(x0-5:x0-1,n))+total(outline(x0+nx:x0+nx+4)))/10.0 


endfor 
sbackground average above box 
for n=y0+ny,139 do begin 
bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 
endfor 
endif else begin 
;box above narcissus 
sbackground average below box 
for n=0,y0-1 do begin 
bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 
endfor 
sbackground average around box 
for n=y0,y0+ny-1 do begin 
bkgdave(n)=(total(outline(x0-5:x0-1,n))+total(outline(x0+nx:x0+nx+4)))/10.0 


endfor 
sbackground average above box 
for n=y0+ny,139 do begin 
bkgdave(n)=(total(outline(0:139,n)))/140.0 
endfor 
endelse 


0 ARES KK a OK RG IK A A A A A A AE KK EE RE EE KE RR REE EBB KEK KE RR EE 
’ 


sgenerate outline and ship average values 
LEAR L ASSES ETE TSR Te TA Te ee ee ene eee ROO eee 
outline(*,0:y0-1)=0.0 
swaterline 
temp1=x0+(9.0/94.0)*nx 
temp2=x0+(89.0/94.0)*nx 
for m=temp1,temp2 do begin 
if outline(m,y0) le BF*bkgdave(y0) then outline(m,y0)=0.0 
endfor 
boat(temp1:temp2,y0)=255 
outline(0:temp1-1,y0)=0.0 
outline(temp2+1:139,y0)=0.0 
;angled portions of bow and stern 








for n=y0+1,y0+(9.0/23.0)*ny do begin 
temp1=x0+(9.0/94.0)*nx-deltal 
temp2=x0+(89.0/94.0)*nx+delta2/3.0 
for m=temp1,temp2 do begin 
if outline(m,n) le BF*bkgdave(n) then outline(m,n)=0.0 
endfor 
boat(x0+(9.0/94.0)*nx-deltal,n)=255 
boat(x0+(89.0/94.0)*nx+delta2/3.0,n)=255 
outline(0:x0+(9.0/94.0)*nx-delta1-1,n)=0.0 
outline(x0+(89.0/94.0)*nx+delta2/3.0+1:139,n)=0.0 
deltal=delta1+1 
delta2=delta2+1 
endfor 
;top of angled portion of bow 
for n=y0+(9.0/23.0)*ny+1,y0+(15.0/23.0)*ny do begin 
temp1=x0+(75.0/94.0)*nx 
temp2=x0+(89.0/94.0)*nx+delta2/3.0 
temp3=x0+(22.0/94.0)*nx 
temp4=x0+(58.0/94.0)*nx 
for m=temp1,temp2 do begin 
if outline(m,n) le BF*bkgdave(n) then outline(m,n)=0.0 
endfor 
for m=temp3,temp4 do begin 
if outline(m,n) le BF*bkgdave(n) then outline(m,n)=0.0 
endfor 
boat(x0+(89.0/94.0)*nx+delta2/3.0,n)=255 
boat(temp3,n)=255 
boat(temp4,n)=255 
boat(temp1,n)=255 
outline(x0+(89.0/94.0)*nx+delta2/3.0+1:139,n)=0.0 
outline(x0+(58.0/94.0)*nx+1:x0+(75.0/94.0)*nx-1,n)=0.0 
outline(0:x0+(22.0/94.0)*nx-1,n)=0.0 
delta2=delta2+1 
endfor 
boat(x0+(75.0/94.0)*nx:x0+(89.0/94.0)*nx+delta2/3.0,y0+(15.0/23.0)*ny)=255 
fantail and foredeck 
temp5=y0+(9.0/23.0)*ny 
boat(x0:x0+(22.0/94.0)*nx,temp5)=255 
boat(temp4:temp1,temp5)=255 
sback of superstructure 
temp1=x0+(22.0/94.0)*nx 
temp2=x0+(58.0/94.0)*nx 
temp4=y0+(15.0/23.0)*ny+1 
temp5=y0+(19.0/23.0)*ny 
temp6=y0+(9.0/23.0)*ny 
for n=temp4,temp5 do begin 
for m=temp1,temp2 do begin 
if outline(m,n) le duck*bkgdave(n) then outline(m,n)=0.0 
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endfor 
endfor 
boat(temp1,temp4:temp5)=255 
outline(0:x0+(22.0/94.0)*nx-1,temp4:139)=0.0 
outline(x0+(58.0/94.0)*nx+1:139,temp4:139)=0.0 
;top of superstructure rear of stack 
temp1=x0+(22.0/94.0)*nx 
temp2=x0+(58.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny 
boat(temp1:x0+(31.0/94.0)*nx,temp3)=255 
srear of stack 
temp1=x0+(31.0/94.0)*nx 
temp2=x0+(44.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny+1 
temp4=y0+ny-1 
for n=temp3,temp4 do begin 
for m=temp1,temp2 do begin 
if outline(m,n) le duck*bkgdave(n) then outline(m,n)=0.0 
endfor 
endfor 
boat(temp1,temp3:temp4)=255 
outline(0:temp1-1,temp3+1:139)=0.0 
outline(temp2+1:139,temp3+1:139)=0.0 
stop of stack 
boat(temp1:temp2,temp4)=255 
outline(*,temp4+1:139)=0.0 
:front of stack: stack average already calculated 
temp1=x0+(44.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny+1 
temp4=y0+ny-1 
boat(temp1,temp3:temp4)=255 
;top of superstructure front of stack: already calculated 
boat(x0+(44.0/94.0)*nx:x0+(58.0/94.0)*nx,y0+(19.0/23.0)*ny)=255 
;front of superstructure: superstructure already calculated 
boat(x0+(58.0/94.0)*nx,y0+(9.0/23.0)*ny:y0+(19.0/23.0)*ny)=255 


a RA OH OR ee ae oF KO OK OK OK OR OB EE EE ROK EK BRE RE REE REE BER EK EK 
’ 


sreposition outline if needed 

AT ORES PERRO Tete eRe T Ter Aer eer eter eer ere 
window,20,xsize=140,ysize=140,xpos=150,ypos=0 

tvscl, boat 

read,'reposition outline? yes=1 no=0',repos 

if repos eq 1 then goto,redo 


oR A HE Ro EK He of a EK a KC A OK OK A EE OK EE OK ORE RE RE EK OR KE ER EE RK 
, 


0 RG A 2H he 2 2 a oO a A A EA EE EE EE KE EER KARE EEE EK KE RK 
? 








scalculate results and display 
EERE EREELECEEER EELS CE TESE ELAR ESE EES CREE ES EEEEEEE 
if image eq 1 then begin 
;save image and average arrays 
bkgdavel=bkgdave 
outlinel=outline 
imagel=boat 
sincrement and go to next image 
image=image+1 
goto,nextimg 
endif 
if image eq 2 then begin 
;save image and average arrays 
outline2=outline 
for n=y0,y0+ny do begin 
shptot1(n)=total(outline1(*,n)) 
shptot2(n)=total(outline2(*,n)) 
index1=where(outlinel(*,n),count) 
shppix1(n)=count 
index2=where(outline2(*,n),count) 
shppix2(n)=count 
if shppix2(n) ne 0.0 then begin 
shpavel(n)=shptot1(n)/shppix2(n) ;use horizontal ship area to find average 
shpave2(n)=shptot2(n)/shppix2(n) 
endif else begin 
shpavel(n)=0.0 
shpave2(n)=0.0 
endelse 
endfor 
bkgdave2=bkgdave 
image2=boat 
;ealculate average ship and bkgd values above and below horizon 
maxi1=max(where(shpavel gt bkgdavel)) 
minl=min(where(shpavel gt bkgdavel1)) 
shpbel1=total(shpavel(min1:yv1))/(yv1-min1+1) 
shpabv1=total(shpavel(yv1+1:max1))(max1-(yv1+1)+1) 
bkgbel1=total(bkgdavel(min1:yv1))/(yv1-min1) 
bkgabv1=total(bkgdavel(yv1+1:max1))/(max1-(yv1+1)+1) 
;display image with outline 
title3=titlel 
window,3,xsize=140,ysize=140,xpos=150,ypos=0,title=title3 
tvscl,imagel 
sgenerate and display outlined image with excess background 
outline1=outline(x0-10:x0+nx+20-1,y0-10:y0+ny+20-1) 
shpimg1=boat(x0-10:x0+nx+20-1,y0-10:y0+ny+20-1) 
title5=title1 
window,5,xsize=nx+20,ysize=ny+20,xpos=300,ypos=0,title=title5 
tvscl,shpimg1 
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;display ship and background average profiles 
title7='ship1' 
window,7,xsize=280,ysize=280 ,xpos=450, ypos=0,title=title7 
plot,shpavel 
oplot,bkgdavel 

scalculate and display contrast profile 
entrstl=(shpavel-bkgdavel)/bkgdavel 
title9='cntrstl' 
window,9,xsize=280,ysize=280,xpos=0,ypos=300,title=title9 
plot,cntrst1(min1:max1) 

scalculate contrast above and below horizon for first image 
entlbel=(shpbell-bkgbel1)/bkgbel1 
entlabv=(shpabv1-bkgabv1)/bkgabv1 

scalculate average ship and bkgd values above and below horizon 
min2=min(where(shpave2 gt bkgdave2)) 
max2=max(where(shpave2 gt bkgdave2)) 
shpbel2=total(shpave2(min2:yv2))/(yv2-min2+1) 
shpabv2=total(shpave2(yv2+1:max2))/(max2-(yv2+1)+1) 
bkgbel2=total(bkgdave2(min2:yv2))/(yv2-min2+1) 
bkgabv2=total(bkgdave2(yv2+1:max2))/(max2-(yv2+1)+1) 

;display image with outline 
title4=title2 
window,4,xsize=140,ysize=140,xpos=150,ypos=160,title=title4 
tvscl,image2 

generate and display outlined image with excess background 
title6=title2 wu 
outline2=outline(x0-10:x0+nx+20-1,y0-10:y0+ny+20-1) 
shpimg2=boat(x0-10:x0+nx+20-1,y0-10:y0+ny+20-1) 
window,6,xsize=nx+20,ysize=ny+20,xpos=300,ypos=160,title=title6 
tvscl,outline2 

generate and display unpolarized outlined image with excess background 
shpimgup=shpimg1/2.0+shpimg2/2.0 
outlinup=outline1/2.0+outline2/2.0 
titleO='unp' 
window,0,xsize=nx+20,ysize=ny+20,xpos=300,ypos=80,title=title0 
tvscl,shpimgup 

display ship and background average profiles 
title8='ship2' 
window,8,xsize=280,ysize=280,xpos=730,ypos=0,title=title8 
plot,shpave2 
oplot, bkgdave2 

calculate and display horizontal contrast profile 
entrst2=(shpave2-bkgdave2)/bkgdave2 
title10='cntrst2’ 
window, 10,xsize=280,ysize=280,xpos=0,ypos=300,title=title10 
plot,entrst2(min2:max2) 

scalculate and display unpolarized contrast profile 

shpaveup=shpavel/2.0+shpave2/2.0 


89 





bkgaveup=bkgdavel/2.0+bkgdave2/2.0 
entrstup=(shpaveup-bkgaveup)/bkgaveup 
title12='entrstup' 

window, 12,xsize=280,ysize=280,xpos=0,ypos=0,title=title12 
plot,entrstup(y0:y0+ny-1) 

;calculate contrast above and below horizon for second image 
cent2bel=(shpbel2-bkgbel2)/bkgbel2 
ent2abv=(shpabv2-bkgabv2)/bkgabv2 

;calculate contrast above and below horizon for unpolarized image 
shpupbel=shpbel1/2.0+shpbel2/2.0 
shpupabv=shpabv1/2.0+shpabv2/2.0 
bkgupbel=bkgbel1/2.0+bkgbel2/2.0 
bkgupabv=bkgabv1/2.0+bkgabv2/2.0 
centupbel=(shpupbel-bkgupbel)/bkgupbel 
entupabv=(shpupabv-bkgupabv)/bkgupabv 

scalculate contrast improvement above and below horizon for v and h 
impbell=cnt1bel/cntupbel 
impabvl=cntlabv/cntupabv 
impbel2=cnt2bel/cntupbel 
impabv2=cnt2abv/cntupabv 

sprint contrast values above and below horizon 
print,'contrast above horizon’ 





print,’ vert hor unp' 
print,cntlabv,cnt2abv,cntupabv 
prin oe 

print,'contrast below horizon 
print,’ vert hor unp' 
print,cnt1lbel,cnt2bel,cntupbel 
prin ie 


print,'contrast improvement above/below horizon 
print,'vertical: ‘',impabvl,impbell 
print, 'horizontal:',impabv2,impbel2 

endif 


sendfor 

goto, loop 
problem: close,1 
print,'I had a problem and so I quit’ 
end 











B. CONTIGUOUS AREA METHOD PROGRAM 


flirashp finds the ship in a picture completely automatic 
scursor is asked for if at picture boundary 

;stores output in chosen file name 

suses shipfind 

‘modular version uses readflir and shipfind 


iscale=2 scale displayed pictures by 2x 
close,1 sclose in case of crash 
on_ioerror,problem sif error, quit program 

path1=' ; 

file1=' ; 


read,'input path to * .img files',path1 

if path1 ne " then path=path1 

print,'path= ',path 

read,'input type of files eg *.img or names to search for',filel 


if filel eq " then filel='*.img 

files=findfile(path+'\'+file1 ,count=num) 

for i=0,num-1 do print,files(i),i 

pene] OgD point to start program over after successful processing**** 


loop: print,n_elements(files)-1,' last number' 

read,'input i for number to plot or -1 to stop’,i 

istart=i 

if i It 0 then goto, problem 

read,'input stop number',istop 

iflfix((istop-istart+1)/2)*2 ne fix(istop-istart+1)) then begin 
print,istart,istop,' must use even number of points' 


goto, loop 
endif 
read,'input first guess difference t level, eg 5 to include all of ship',tdif 
tdifsave=tdif ssave tdif so that changes will not affect batch 


ie aac ar aia sada (016) POUL SPST eer E Te eae Sees 
contrast=dblarr((istop-istart+1)/2,5) 
for ii=istart, istop,2 do begin 

for i=ii,ii+1 do begin 


tdif=tdifsave -restore the tdif if changed in loop 

f=files(i) 

file=f 

readflir,f,yy,t,bar,tt,head smodular input returns binary, temparature 


SEMEL EEE TET PER ER SSeS display data in windows for visual inspection *********** 


forder=0 











ct=0 

loadct,ct 

if(ii eq i) then window,0,xsize=280,ysize=280,xpos=0,ypos=0,title='0 binary data’ $ 
else window,5,xsize=280,ysize=280,xpos=280,ypos=0,title='5 binary unscaled data’ 


print,'UNScaled plot of binary data in 0' 
;tvscl,congrid(yy,140*iscale,140*iscale,Interp=1) 
tv,congrid(yy,140*iscale,140*iscale,Interp=1) 


wait,1 
loadct,16 
window,1,xsize=280,ysize=330,xpos=600,ypos=0,title='1 temperature’ 


tvscl,rebin(tt,140*iscale,165*iscale) 


eke em perature point option bypassed ****## #2 i a He He I ee ea 


goto,skip 

print,'Click in scale bar to exit’ 

print,'click in picture for point temperature’ 
replot: cursor,xv,yv,3,/device 

if yv lt 50 then goto,skip 
print,t(xv,330-yv),mint,maxt 

goto, replot 

skip: 


0 FCG ae 2 2 ee he oe oe oe kok ae ok i AK a ae ok ACK Ak AK AE oO] Shai FER RE aE eR eb eo 
’ 


RARER ERATE EEATEE AO two bars at right and left of picture and 


REREEEEEAEEEREREREREREE EY SO the lower average value, then smooth it 


thresO=(t(137,*)+t(138,*)+t(139,*))/3.0 
thres=(t(0,*)+t(1,*)+t(2,*))/3.0 
thres=thres > thres0 ;changed to larger of backgrounds 


minthres=min(thres) 

maxthres=max(thres) 

avthres=total(thres)/n_elements(thres) 
thres=smooth(thres,3) 


thres=rebin(thres, 140,140) smake the vertical line xy array 
thresO=thres ;save in thresl for debug 140x140 array 
sship=(t gt thres+tdif)*t sship is greater than threshold 


sshipb=(t gt thres+tdif}*yy get binary version for display 





ship=(t gt thres)*t ;ship is greater than threshold 


shipb=(t gt thres)*yy sget binary version for display 
snew version uses threshold tdif later 


if(total(ship) eq 0) then stop ;zero values mean error< 


;***make linear array for averages 
shipt=ship(where(ship gt 0)) _ ;Pick only nonzero values for average 


;sshipt is linear array now 
avship=total(shipt)/n_elements(shipt) 
maxship=max(shipt) 
minship=min(shipt) 
-print,minship,maxship,avship,' ship min,max,average temperature’ 


:rint,minthres,maxthres,avthres,' min max average background! 


‘Find the MAXIMUM temperature 


Point ee renner Teen een nee 


maxt=max(ship,it) ;find the max in the ship thresholded 


ix=it mod 140 
iy=it/140 


window,4,xsize=280,ysize=280,xpos=600,title='#4 binary ship' 
tvscl,rebin(shipb,140*iscale,140*iscale) smake ship twice as big for display 


tvers,ix*iscale,iy*iscale,/device sput the cursor on the peak location 


:print,'right button, use auto position, left, new cursor position’ 

lerr=4 

if(i ne istart) then begin 
print,ix,iy,nx,ny,x0,y0,format='(6(i3,2x),A19)',' ix,iy,nx,ny,x0,y0' 


iffix 1t 10 or iy It 10) then cursor,x0,y0,/device sAUTOMATIC LINE 
OUT 
if(nx/ny gt 5 or ny It 10 or nx lt 10) then cursor,x0,y0,/device 
‘AUTOMATIC LINE OUT ;left button 1, middle 2 right 4 
endif 
if(lerr eq 1) then begin 
iy=y0/iscale 
ix=x0O/iscale 
xO=ix.. 








yO=iy 


endif else begin 
x0=1x 
yO=iy 


endelse 
SEES EEERERE EERE ECE ETERS ID loop find POINT Fe TRS Ser eS See Tee TEER Re reer 


:following necessary for reship loop, otherwise ship will not change with threshold 


reship: 
ESE eee net hind 
C1 lle 


x0=ix 

yO=iy 
shipfind,shipb,ship,x0,y0,nx,ny,iscale use find on the binary version to avoid neg t 
;WARNING SHIPFIND SETS ALL VALUES IN SHIP AND SHIPB ZERO NOT ON 


THE SHIP 
tvscl,rebin(shipb,140*iscale,140*iscale) ;show the contiguous thresholded=0 ship 


sbox:xy=boxfun(x0,y0,nx,ny,/init,/message);option to readjust the box 


boxdraw,x0,y0,nx,ny ;show where the box is 
x0=x0/iscale srescale from picture to data size 
y0=y0/iscale 


nx=nx/iscale 
ny=ny/iscale 


;print,x0,y0,nx,ny 

thres=(t(x0+nx,*)) ‘right side of box 

thres1=t(x0,*) sleft side of box 

thres=thres > thres1 ;changed to greater than for background 


thres=smooth(thres,2) 

thres=rebin(thres, 140,140) smake an array 140x140 
thres1=fltarr(140,140) 

thres1(*,*)=0.0 

thres1(x0:x0+nx,y0:y0+ny)=thres(x0:x0+nx,y0:y0+ny) ;extract part of array 
;Now use the threshold to clean up the ship and check for loop 

:the ship and shipb cannot be any larger than tdif=0, so whittle 

ship=(t gt thres1+tdif)*ship ;ship is greater than threshold 


shipb=(t gt thres1+tdif)*shipb 


;get binary version for display 
xO=ix 
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yO=iy 
shipfind,shipb,ship,x0,y0,nx,ny,1 suse find on the binary version to avoid neg t 


window,2,xsize=280,ysize=280,xpos=300,title='#2 local threshold’ 
tvscl,rebin(shipb,140*iscale,140*iscale) sredisplay 
sabove is entire 140x140 picture thresholded 
box=t(x0:x0+nx,y0:y0+ny) ;temperature in box 
raw=yy(x0:x0+nx,y0:y0+ny) sraw signal in box 
ship1=ship(x0:x0+nx,y0:yO+ny) ;ship temperature in box 
raw1=shipb(x0:x0+nx,y0:yO+ny) ;ship binary in box 
shiplb= box*(ship1 eq 0) ;ship temperature BACKGROUND in box 


rawlb= raw*(raw1 eq 0) sship binary BACKGROUND in box 
;ship1 and raw1 are box sized subset of ship and shipb 


window,3,xsize=nx*iscale*2, ysize=ny*iscale*2,xpos=0,ypos=0 
a> d > 


tvscl,congrid(ship1,(nx+1)*iscale*2,(ny+1)*iscale*2) 


tdif1=0.0 
sread,'input new tdif or 0 to stop ',tdifl ;AUTOMATIC LINE OUT 
if(tdifl gt 0) then begin 
tdif=tdifl 
remake the ship variables using the NEW box threshold 
ship=(t gt thres+tdif)*t ;ship is greater than threshold 
shipb=(t gt thres+tdif)*yy ;get binary version for display 
goto,reship sloop up to reship*#*## +++ # HEE HARE AREA AAA AA ALAR 
endif 


POINT 88 88 F EE EERE REAR RH EE A BR IR HR TER ER EE EE ET REE 
REET ETEEER TESTES PLONE RES CATES ERE TSave 
SCIP ARE BBE EERE oii sok iio eR RSI ICR CE 
if(i eq ii) then begin 
tdifsave=tdif 


xOsave=x0 ;box position 
yOsave=y0 
nxsave=nx 
nysave=ny 
ixsave=ix speak position 
iysave=ly speak position 


endif 


ifGi eq ii) then begin 








shiplsv=ship1 
rawlsv=rawl 
shiplbsv=ship1b 
raw lbsv=raw1b 


tsave=t 

yysave=yy 

shipsv=ship ssave the temperature thresholded full picture 
shipbsv=shipb ;save the binary thresholded 


endif else begin 


sif the first box is bigger, renormalize to it 
;otherwiser renormalize all data to the second box 
;The following are necessary for correct contrast using binary 
sraw1,rawlb,rawlsv,rawlbsv - only two of these are destroyed 
sonly one set of the two needs to be recomputed based on the new box 
;size which is larger 
print, nxsave,nysave,nx,ny,' nxsave,y,nx,ny’ 
if(nxsave*nysave gt nx*ny) then begin 
print,'using old ship size’ 
‘use save box and make raw1 and rawlb match old ship shipsv 
x0=x0save 
y0=y0save 
nx=nxsave 
ny=nysave 


rawl=yy*(shipbsv gt 0) ;mask saved correct ship on current data 
rawlb=yy*(shipbsv eq 0) smask for background = not ship 
raw1=raw1(x0:x0+nx,y0:y0+ny) sselect the current (saved ) box 
raw1b=raw1b(x0:x0+nx,y0:y0+ny) / 


endif else begin 
print,'using new ship size’ 
rawlsv= yysave*(shipb gt 0) smask current ship on old data 140x140 
raw lbsv=yysave*(shipb eq 0) 
;current not ship binary times old data 
raw 1sv=raw1sv(x0:x0+nx,y0:y0+ny) j;now select the current box 
raw lbsv=raw1bsv(x0:x0+nx,y0:yO+ny) ;background 
sit would also be nice later to update the ship1 (temperature) stuff 
endelse 
endelse 


aOR IE A A A OE EK KK OK AK KE IK OK A EK KE EO IR RR EK BE AE EE RE EE 
) 


endfor 
eC CEE A MH ER KE OK RR EE RE EE RE HK ER EERE EE RK EEE EK EE KK EK 
? 


RK 


96 





EPELERE AEDT ECES RESTA EEE ERET ESE icnlay* ttc ee ee ee eee 
, 

ORE 

able ask a oleae 2 of eae a oR KC aR OR a AC a AACE OE a a EH HE AA OR Ho OI HE OR EH HA RAE EA HE HE EE AA a IC OR OR SEH HO OE A OK KO KR OR 
; i 

SOR 


loadct,0 

sin order 

sship 1 average value 

;sship 1 background average value 
;ship 2 average value 

;ship 2 background average value 
sship unpolarized 

ship background unpolarized 

sship unpolarized average value 
sship background unpolarized average value 
scontrast shipl 

contrast shipsv (ship0) 

scontrast unpolarized 

scontrast shipi/contrast unpolarized 
;contrast shipsv/contrast unpolarized 


shiplav=total(ship1)/total((ship1 gt 0)) 

shiplbav = total(ship1b)/total((ship1b gt 0)) 
shipsav=total(ship1sv)/total( (ship1sv gt 0)) 

shipsbav = total(ship1lbsv)/total( (ship1bsv gt 0)) 
shipup=(ship1+ship1sv)/2.0 

shipbup=(ship1b+ship1bsv)/2.0 
shipupav=total(shipup)*2.0/total((ship1 gt 0)+(shiplsv gt 0)) 
shipubav=total(shipbup)*2.0/total((ship1b gt 0)+(ship1sv gt 0)) 


‘above are taken for temperature, below for raw dataq 


shiplav=total(raw1)/total((raw1 gt 0)) 

shiplbav = total(raw1b)/total((raw1b gt 0)) 
shipsav=total(raw1sv)/total( (raw1sv gt 0)) 

shipsbav = total(rawlbsv)/total( (rawlbsv gt 0)) 

;WARNING, adding binary will exceed 255 and give error 

sone must first float the numbers; using temperature above is ok 
shipup=(float(raw1)+float(raw1sv))/2.0 
shipbup=(float(raw1b)+float(raw1bsv))/2.0 
shipupav=total(shipup)*2.0/total((raw1 gt 0)+(rawlsv gt 0)) 
shipubav=total(shipbup)*2.0/total((raw1b gt 0)+(rawlbsv gt 0)) 


;total shipup gt 0 does not work because of registration 


eship1=(shiplav-shiplbav)/shiplbav 
cshipsv=(shipsav-shipsbav)/shipsbav 
eshipu =(shipupav-shipubav)/shipubav 
cclu=cship1/cshipu 











ecsu=cshipsvw/cshipu 

srawl = ship picture in binary 0-255 (in box) 
srawlsv= same picture 0 

srawlb = binary background in box 
srawlbsv = binary saved background 

syy = full picture (same as t for temperature 
syysave= picture 0 binary 

;t=temperature picture 


;tsave=old picture temperature 
0 RO OK ee he oie of ae KK OR a KR EK a EK OK KK ER RE KE RRA REE E REE RRR E REE EK 
a 


Print,'Contrast and enhancement,first file =',file,tdif,tdifsave 


iftistop-istart gt 0) then begin 
print,files(ii+1),files(i),' unpolarized ,first/unpolarized,second/unpolarized' 


print,cship1,cshipsv,cshipu,cclu,ccsu 
contrast((ii-istart)/2,0:4)=[cship1,cshipsv,cshipu,cclu,ccsu] 
ssave the values 


endif else begin 

sprint,’ first contrast ',file 
;print,cship1 

endelse 


display difference image for both polarizations 
window,2,xsize=nx*4,ysize=ny*4 xpos=0,ypos=0,title=files(istop) 


tvscl,rebin(raw1,(nx+1)*4,(ny+1)*4) 
window,3,xsize=nx*4,ysize=ny*4,xpos=0,ypos=140,title=files(istart) 


tvsel,congrid(rawl1sv,(nx+1)*4,(ny+1)*4) 
AR ELEASE CET ERE EEE EE EE TEE TEE ERTS EPS Sei plots in automatic 
MOE TEE EEE AREER AAA AER RR AR RI 
goto, noplot 
:plot ship average value compared to background average value and make .tif file 


shorizontal polarization 
titlel="horizontal'+' polarization’ 
window,4,xsize=350,ysize=300,xpos=300,ypos=0,title=titlel __ ;files(istop) 


wset,4 
scall to external function avrow averages rows 
plot,avrow(ship1sv),indgen(140),background=255,color=0 


oplot,avrow(ship1bsv),indgen(140),color=150 


; oplot,bkgdave(*),indgen(140),color=150 
ploti=tvrd() 
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plot1=plot1(*,299-indgen(300)) 
tiff_write,title1+'.tif,plot1 


splot comparison of ship average values with respect to polarization and make .tif file 


title2='polarization'+' comparison’ 
window,5,xsize=350,ysize=300,xpos=300,ypos=330,title=title2 ;'vertical/horizontal' 


wset,5 
plot,avrow(ship1),indgen(140),background=255,color=0 


oplot,avrow(ship1sv),indgen(140),color=150 
plot2=tvrd() 
plot2=plot2(*,299-indgen(300)) 
tiff_write,title2+'.tif ,plot2 


:plot first ship average value compared to background average value and make .tif file 


svertical polarization 
title3='vertical'+' polarization’ 
window, 6,xsize=350,ysize=300,xpos=660,ypos=0,title=title3 _;files(istart) 


wset,6 
plot,avrow(ship1),indgen(140),background=255,color=0 


oplot,avrow(ship1b),indgen(140),color=150 
plot3=tvrd() 
plot3=plot3(*,299-indgen(300)) 
tiff_write,title3+'.tif ,plot3 


splot comparison of backgrounds with respect to polarization and make tif file 


title4='background'+' comparison’ 
window,7,xsize=350,ysize=300,xpos=660,ypos=330,title=title4 


wset,7 
plot,avrow(ship1bsv),indgen(140),background=255,color=0 


oplot,avrow(ship1b),indgen(140),color=150 
plot4=tvrd() 
plot4=plot4(*,299-indgen(300)) 
tiff_write,title4+'.tif plot4 
URE ET EEE RT Ce Tee ee end of temperature 
WIN COWS EH ERER ER ER ERA EAA A EA RARE EE ERE ERE 


:plot ship average value compared to background average value and make .tif file 


shorizontal polarization 
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title1='hzpol' 
window, 10,xsize=350,ysize=300,xpos=300,ypos=0,title=titlel _;files(istop) 





seall to external function avrow averages rows 
plot,avrow(raw1sv),indgen(140),background=255,color=0,title="horizontal'+' polarization’ 


oplot,avrow(raw1bsv),indgen(140),color=150 
; oplot,bkgdave(*),indgen(140),color=150 
plotl=tvrd() 

plot1=plot1(*,299-indgen(300)) 
tiff_write,title1+'.tif,plot1 


splot comparison of ship average values with respect to polarization and make .tif file 


title2='comp' 

window, 12,xsize=350,ysize=300,xpos=300, ypos=330,title=title2 
;'vertical/horizontal' 

plot,avrow(raw1),indgen(140),background=255,color=0,title='polarization'+' 
comparison' 


oplot,avrow(raw1sv),indgen(140),color=150 
plot2=tvrd() 
plot2=plot2(*,299-indgen(300)) 
tiff_write,title2+'.tif ,plot2 


splot first ship average value compared to background average value and make .tif file 


;vertical polarization 
title3='vpol' 
window, 13,xsize=350,ysize=300,xpos=660,ypos=0,title=title3 __;files(istart) 


plot,avrow(raw1),indgen(140),background=255,color=0,title='vertical'+' polarization’ 


oplot,avrow(raw1b),indgen(140),color=150 
plot3=tvrd() 
plot3=plot3(*,299-indgen(300)) 
tiff_write,title3+'.tif plot3 


splot comparison of backgrounds with respect to polarization and make .tif file 


title4='bkgnds' 
window, 14,xsize=350,ysize=300,xpos=660, ypos=330,title=title4 


plot,avrow(raw1bsv),indgen(140),background=255,color=0,title="background'+' 
comparison’ 

oplot,avrow(raw1b),indgen(140),color=150 

plot4=tvrd() 
plot4=plot4(*,299-ir. an(300)) 











tiff_write,title4+'.tif ,plot4 


noplot: 
endfor 


;saveflr8.pro prints results 

;saveflr8.pro prints results 

fileout=' ; 

close,2 

read,'input output filename’ fileout 

openw,2,fileout 

printf,2,systime(),' date of analysis’ 

printf,2,'Temperature difference = ',tdif 
printf,2,istop-istart+1,'number of files' 

j=istart 

for i=0,(n_elements(contrast(*,0))-1) do begin;printf,2,files(2*i),files(2*i+1) 
printf, 2, files(j+1), files(j),format='(2(A30))’ 

j=j+2 

printf,2,'first,second,unpolarized ,first/unpolarized,second/unpolarized' 
printf,2,contrast(i,*),format="(5(f9.6,5x))" 
print,i,contrast(i,*),format="(i3,1x5(f9.6,5x))" 

endfor 


close,2 
print, 'finished’ 


goto, loop 

problem: close,1 

print,'I had a problem and so I quit’ 
end 


C. BOX METHOD PROGRAM 


slasthope.pro 
lorder=0 
print,'flirs makes movie of flirs' close,1 
sclose in case of crash intfn=0 
shpnmbr=1 
ct=0 loadct,ct 
;window,0,xsize=280, ysize=280,xpos=0,ypos=0 
on_ioerror,problem path1=' 
, filel=' 
' print,'file size must be 20446 
with header' read,'input path to *.img files',path1 


if path1 ne " then path=path1 
print,'path= ',path 
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read,'input type of files eg *.img or names to search for',filel 
if filel eq " then filel='*.img 
files=findfile(path+' \'+file1,count=num) for 
i=0,num-1 do print,files(i),i 


0 DESI HER OR A AK oH A RK A KK EK A KKK EK AR OR EE RK EEE EEK REE RE RE RK 
, 


sprint time of analysis and number of files analyzed 

PELLET TTS ete ee SORT ENE ee CELE E TELA R ENE Se ST ETE CE ee hee 
printf,7,' ‘ systime() 

printf,7,' ' strmid(files(5),17,5),'.txt' 


0 RS C2 a ae 2 a 2 oS A A A RK OK A EK KE EE KE AER RE ER OK EEK RE KEE KE EE RE 
? 


loop: print,n_elements(files)-1,' last number' 


read,'input i for number to plot or -1 to stop',i ifilt 0 
then goto, problem istart=i 
read,'input stop 
number',istop for i=istart, istop do begin 
f=files(i) 


ORR RRR EEE RAKE REG RRR RE RA RRREECEE E EE 
’ 


;open output file 

EET ETAT UREA S ESTATE RESL STP RET ERT STC ee Nene te ene eee 
openr,1,f 

print,'opening ',f,i 

head=bytarr(846) 

y=bytarr(140,70) 

yl=bytarr(140,70) 

yy=bytarr(140,140) 


readu,1,head readu,l,y 
readu,l,yl 


close,1 
index=indgen(70) 
yy(*,index*2)=y1(*,index) 
yy(* ,index*2+1)=y(*,index) 
if(max(yy) eq 255) then begin 
yl=yy( where(yy gt 254)) 


nover=n_elements(y1) endif else begin 
nover=0 


endelse 
yy=yy(*,139-indgen(140)) 


0 RSI KK A EE EK EE KE KEE REE RE KK KK 
’ 


;save ship images 

EERE R ECON T LR nae CCT Tee 
shipl=bytarr(140,140) 
ship2=bytarr(140,140) 

if shpnmbr eq 1 then begin 
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shipl=yy 
shipll=yy 
shpnmbr=2 
shpname1=files(i) 
titlel='vert' 
window, 1,xpos=0,ypos=0,xsize=140,ysize=140,title=title1 
tvscl,ship1 
endif else begin 
ship2=yy 
shpname2=files(i) 
title2="hor' 
window,2,xpos=0,ypos=160,xsize=140, ysize=140, title=title2 
tvscl,ship2 
shpnmbr=1 
endelse 
endfor 


See nvert sional stFeernr yl=1 
y=1 


KEE CEE pO cot variables*****####* 


yy=congrid(yy,140*2,140*2,Interp=1) ;lorder=1 


loadct,0 
swset,0 
print,'Scaled plot of data in 0' stvscl yy 
range=double(head,770) 

level=double(head,762) 

dlevel=double(head,618) 

emiss=double(head,626) 

ct=double(head,650) 

r=double(head,570) 

b=double(head,578) 

f=double(head,586) 

amt=double(head,674) 

at=double(head,666) 

pl=emiss*ct 

p2=(1-emiss)/emiss 

p3=1.0-ct 

iamb=r/(exp(b/amt)-f) 

iatm=r/(exp(b/at)-f) 

y=range*(float(yy)-128.)/256.0+level+dlevel 

y=y/p1-p2*iamb-p3/p1*iatm 

t=b/alog(abs((r/y+f)))-273.15 

idate=head(482:487) 
itime=head(488:495) 

; get rid of old variables 
sdelvar,y,yl,yy,head 








shpup=fltarr(140) 
shpup(*)=0.0 
shp1=fltarr(140) 
shp1(*)=0.0 
shp2=fltarr(140) 
shp2(*)=0.0 
bkgdup=fitarr(140) 
bkgdup(*)=0.0 
bkgd1=fltarr(140) 
bkgd1(*)=0.0 
bkgd2=fltarr(140) 
bkgd2(*)=0.0 
shpave=fltarr(140) 
shpave(*)=0.0 
sshpavel=fitarr(140) 
;shpavel(*)=0.0 
sshpave2=fltarr(140) 
;shpave2(*)=0.0 
;shptot1=fltarr(140) 
;shptot1(*)=0.0 
sshptot2=fltarr(140) 
;shptot2(*)=0.0 
shppix1=fltarr(140) 
shppix1(*)=0.0 
shppix2=fltarr(140) 
shppix2(*)=0.0 
outline=fitarr(140,140) 
outline(*,*)=0.0 
outlinel=fltarr(140,140) 
outline1(*,*)=0.0 
outline2=fltarr(140,140) 
outline2(*,*)=0.0 


image=1 

nextimg: 

wset,image 
zz=boxfun(x0,y0,nx,ny,init=intfn,/message) 
intfn=1 

box=tvrd(x0,y0,nx,ny) 

bkgdave=fltarr(140) 

bkgdave(*)=0.0 

shppix=fltarr(140) 

shppix(*)=0.0 


eI AE RE ER EK KE EE REE ERK REE EKER ERE EERE KK 
? 


assign image array to outline array and mark horizon 
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a 


if image eq 1 then begin 
outline=ship11 
print,'click on horizon or bottom to exit’ 
replotl: cursor,xv,yv,3,/device 
if yv lt 5 then goto,exit1 
xvl=xv 
yvl=yv 
YVY=YV 
print,‘horizon=',yv1 
print,'click on point or bottom to exit’ 
goto,replot1 
exit1:print,'exiting’ 
endif 
if image eq 2 then begin 
outline=ship2 
print,'click on horizon or bottom to exit’ 
replot2: cursor,xv,yv,3,/device 
if yv It 5 then goto,exit2 
XV2=xV 
yv2=yv 
YVY=YV 
print,horizon=',yv2 
print,'click on point or bottom to exit’ 
goto,replot2 
exit2:print,'exiting' 
endif 


0 Be oe 246 2s fe 2S 2S a EE A OK Oe oe 2 He 3 3 2 he a OE AE a OK EO A OE A OK KK RE BR REE RE KE EE 
? 


scalculate background average values 
0H a fo a i 2K OK i 2 i 2 OK A a OK KK EE KK EK RK EEK EKER EER KR KE 
? 


if yO+ny/2.0 lt 69 then begin 
sbox below narcissus 
sbackground average below box 
for n=0,y0-1 do begin 
bkgdave(n)=total(outline(0:139,n))/140.0 
endfor 
sbackground average around box 
for n=y0,y0+ny-1 do begin 
bkgdave(n)=(total(outline(x0-5:x0-1,n))+total(outline(x0+nx:x0+nx+4)))/10.0 


sbackground average above box 
for n=y0+ny,139 do begin 
bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 


endif else begin 
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endfor 


endfor 





;box above narcissus 

sbackground average below box 

for n=0,y0-1 do begin 
bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 endfor 





sbackground average around box 


for n=y0,y0+ny-1 do begin 
bkgdave(n)=(total(outline(x0-5:x0-1,n))+total(outline(x0+nx:x0+nx+4)))/10.0 endfor 


sbackground average above box 

for n=y0+ny,139 do begin 
bkgdave(n)=(total(outline(0:139,n)))/140.0 
endfor 

endelse 


0 DG AG a a ae CO ae 2 OK A OK OK AK EK IK EE EE BR EE ERE RE KEE EK 


read,'enter background factor',dd 
duck=dd 
if image eq 1 then begin 
factorl=duck 
endif else begin 
factor2=duck 
endelse 


again: 
shppix(*)=0 
outline(0:139,0:y0-1)=0.0 
outline(0:139,y0+ny:139)=0.0 
outline(0:x0-1,0:139)=0.0 
outline(x0+nx:139,0:139)=0.0 
for n=y0,y0+ny-1 do begin 
for m=x0,x0+nx-1 do begin 
if outline(m,n) le duck*bkgdave(n) then outline(m,n)=0.0 
if outline(m,n) gt 0.0 then shppix(n)=shppix(n)+1 
endfor 
if shppix(n) eq 0 then shpave(n)=0.0 
if shppix(n) ne 0 then shpave(n)=total(outline(*,n))/shppix(n) 
endfor 
numbel=total(shppix(y0:yvy)) 
numabv=total(shppix(yvy+1:y0+ny-1)) 


shpabv=total(outline(x0:x0+nx-1,yvy+1:y0+ny-1)) 
shpbel=total(outline(x0:x0+nx-1,y0:yvy)) 
bkgdabv=total(bkgdave(yvy+1:y0+ny-1))/((yO+ny-1)-(yvy+1)+1) 
bkgdbel=total(bkgdave(y0:yvy))/(yvy-y0+1) 


window,3,xsize=140,ysize=140,xpos=420,ypos=0 
tvscl,outline 
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read,'enter new background factor or zero to continue’,qq 
duck=qq 
if duck gt 0.0 then begin 
if image eq 1 then begin 
outline=ship11 
factorl=duck 
goto,again 
endif else begin 
outline=ship2 
factor2=duck 
goto,again 
endelse 
endif 
if image eq 1 then begin 
outlinel=outline 
shp1(*)=shpave(*) 
bkgd1(*)=bkgdave(*) 
shpabvl=shpabv 
shpbell=shpbel 
bkgdabvl=bkgdabv 
bkgdbell=bkgdbel 
numabvl=numabv 
numbel1l=numbel 
image=image+1 
window,4,xsize=140,ysize=140,xpos=280,ypos=0 
outline(*,yv1)=255 
tvscl,outline 
window,5,xsize=140,ysize=140,xpos=140,ypos=0 
tv,ship11 
goto,nextimg 
endif 
if image eq 2 then begin 
outline2=outline 
shp2(*)=shpave(*) 
bkgd2(*)=bkgdave(*) 
bkgdup(*)=bkgd1(*)/2.0+bkgd2(*)/2.0 
shpup(*)=shp1(*)/2.0+shp2(*)/2.0 


shpabv2=shpabv 
shpbel2=shpbel 
bkgdabv2=bkgdabv 
bkgdbel2=bkgdbel 


shpabvla=shpabv1/numabv1 
shpbella=shpbeli/numbell 


numabv2=numabv 
numbel2=numbel 








if numabv2 ge numabv1 then begin 
shpabv1=(shpabv1+bkgdabv1*(numabv2-numabv1))/numabv2 
shpabv2=shpabv2/numabv2 

endif else begin 
shpabvl=shpabv1/numabv1 
shpabv2=(shpabv2+bkgdabv2*(numabv1-numabv2))/numabv1 

endelse 

if numbel2 ge numbell then begin 
shpbell=(shpbell+bkgdbel1*(numbel2-numbel1))/numbel2 
shpbel2=shpbel2/numbel2 

endif else begin 
shpbell=shpbel1/numbel1 
shpbel2=(shpbel2+bkgdbel2*(numbel1-numbel2))/numbell1 

endelse 


entabv1=(shpabv1-bkgdabv1)/bkgdabv1 
cntbell=(shpbell-bkgdbel1)/bkgdbel1 
cntabv2=(shpabv2-bkgdabv2)/bkgdabv2 
cntbel2=(shpbel2-bkgdbel2)/bkgdbel2 


shpabvu=shpabv1/2.0+shpabv2/2.0 
shpbelu=shpbel1/2.0+shpbel2/2.0 
bkgdabvu=bkgdabv1/2.0+bkgdabv2/2.0 
bkgdbelu=bkgdbel1/2.0+bkgdbel2/2.0 
entabvup=(shpabvu-bkgdabvu)/bkgdabvu 
cntbelup=(shpbelu-bkgdbelu)/bkgdbelu 


impabvh=cntabv2/entabvup 
impbelh=cntbel2/cntbelup 
impabvv=cntabv1/cntabvup 
impbelv=cntbel1/cntbelup 


window,6,xsize=140,ysize=140,xpos=280,ypos=160 
outline(*,yv2)=255 
tvscl,outline 
window,7,xsize=140,ysize=140,xpos=140,ypos=160 
tv,ship2 

endif 


print,' shpbella shpabvla 
print,shpbella,shpabvla 

print,' shpbell shpbel2 shpbelu' 
print,shpbel1,shpbel2,shpbelu 

print,' bkgdbell bkgdbel2  bkgdbelu' 
print, bkgdbel1,bkgdbel2,bkgdbelu 
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print," 

print,’ cntabvl cntabv2 entabvup' 
print,cntabv1,cntabv2,cntabvup 

print,’ cntbell entbel2 cntbelup' 
print,cntbell,cntbel2,cntbelup 

print,’ numabv1l numabv2' 
print,numabv1,numabv2 

print,' numbell numbel2' 
print,numbel1,numbel2 

print,' impabvv impabvh' 
print,impabvv,impabvh 

print,' impbelvy = impbelh' 
print,impbelv,impbelh 

print,'bkgd factor v/h',factor1,factor2 
print,‘horizon v/h',yv1,yv2 


ship=string(‘ship') 
bkgd=string(‘bkgd') 
entabv=string(‘cntabv') 
cntbel=string(‘cntbel’) 
numabv=string(‘numaby') 
numbel=string(‘numbel') 
impabv=string(‘impabv') 
impbel=string(‘impbel') 
BF=string('BF") 
hrzn=string(‘hrzn') 


printf,7,' ' shpname!1,'\',strmid(shpname2,17,12) 

printf,7,'vert',‘hor','unpol' itime(0),itime(2),itime(4),format='(20x,a4,5x,a3,6x,a5,2x,i2,1h:, 
i2,1h:,i2)' 

printf,7,'shipabv' shpabv1,shpabv2,shpabvu,format='(10x,a7,1x,f6.2,3x,f6.2,3x,f6.2)' 
printf,7,'bkgdabv' ,bkgdabv1,bkgdabv2, bkgdabvu,format='(10x,a7,1x,f6.2,3x,f6.2,3x,f6.2)' 


printf,7,'shipbel' shpbel1,shpbel2,shpbelu,format='(10x,a7,1x,f6.2,3x,f6.2,3x,f6.2)' 
printf,7,'bkgdbel' ,bkgdbell1 ,bkgdbel2,bkgdbelu,format='(10x,a7,1x,f6.2,3x,f6.2,3x,f6.2)' 
printf,7,'cntabv' cntabv1,cntabv2,cntabvup,format='(10x,a6,2x,f6.2,3x,f6.2,3x,f6.2)’ 
printf,7,’cntbel' cntbell,cntbel2,cntbelup,format='(10x,a6,2x,f6.2,3x,f6.2,3x,f6.2)' 
printf,7, numabv' numabv1,numabv2,format='(10x,a6,2x,f6.2,3x,f6.2)' 

printf,7,'numbel' numbell numbel2,format='(10x,a6,2x,f6.2,3x,f6.2)' 

printf,7,'impabv' ,impabvv,impabvh,format='(10x,a6,2x,f6.2,3x,f6.2)' 

printf,7 ,'impbel' impbelv,impbelh,format='(10x,a6,2x,f6.2,3x,f6.2)' 

printf,7,'BF" ,factor1,factor2,format='(10x,a2,6x,f6.2,3x,f6.2)' 

printf,7,"hrzn' ,yv1,yv2,format='(10x,a4,4x,f6.2,3x,f6.2)' 

printf,7," 


intfn=0 
shpnmbr=1 








goto, loop problem: 
close,1 print,'T had a 


problem and so I quit’ 
end 
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